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FOREWORD 


This  report  is  one  of  two  computer  program  user's 
manuals  prepared  by  Aerotherm  Corporation  under  USAF 
Contract  F04611-70-C-0012 .  Included  herein  is  Volume 
I  of  the  manual  for  Version  3  of  the  Aerotherm  Equilib¬ 
rium  Surface  Thermochemistry  code.  This  volume  describes 
the  problems  solved  by  the  code  and  presents  an  input 
(card  format)  user's  guide  and  sample  problems.  The 
report  was  first  published  as  Aerotherm  Report  No.  UM-70-13, 
April  1970.  The  work  was  administered  under  the  direction 
of  the  Air  Force  Rocket  Propulsion  Laboratory,  Motor 
Component  Development  Branch  with  Mr.  R.  J.  Schoner  as 
project  officer. 

Mr.  M.  R.  Wool  was  program  manager  and  principal 
investigator.  Significant  additional  assistance  was 
also  provided  by  Dr.  C.  B.  Moyer. 

This  technical  report  has  been  reviewed  and  is 
approved . 


R.  J.  Schoner 
Project  Engineer,  AFRPL 


ABSTRACT 


A  Fortran  IV  computer  code  is  described  which  computes 
the  equilibrium  or  mixed  frozen/equilibrium  thermodynamic 
state  of  general  chemical  systems.  Closed  (fixed  mass)  molec¬ 
ular  compositions  are  evaluated  from  relative  elemental  quan¬ 
tities,  species  thermochemical  data  and  two  state  properties 
(pressure  and  temperature  or  enthalpy) .  Open  system  (diffusive 
mass  flux  dependent)  states  are  defined  from  surface  equilib¬ 
rium  considerations  utilizing  a  film  coefficient  model  which 
accounts  for  the  unequal  diffusion  of  species.  The  generality 
of  the  formulation  allows  computations  for  a  broad  range  of 
environments  and  surface  materials  including  charring  ablators. 
This  computer  code  is  designated  Version  3  of  the  Aerotherm 
Equilibrium  Surface  Thermochemistry  code  (EST3)  and  provides 
surface  mass  balance  quantities  needed  for  ablation  predictions 
by  the  Aerotherm  Charring  Material  Ablation  code  (CMA) . 

Volume  I  of  this  report,  presented  herein,  contains 
descriptions  of  the  fundamental  physical  events  modeled, 
the  mathematical  equations  solved,  the  information  required 
for  input,  and  the  results  output  by  the  computer  code.  An 
input  (card  format)  user's  guide  is  provided  along  with  sample 
input  and  output  listings  to  enable  an  unfamiliar  user  to 
successfully  operate  the  code  and  understand  the  results. 

Volume  II  of  this  report  contains  supplemental  information  on 
the  specific  Fortran  IV  codings.  Included  are  program  listings, 
flow  charts,  and  definitions  of  Fortran  variables. 
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SECTION  1 
INTRODUCTION 


The  computer  program  described  in  this  user's  manual 
is  a  revised  edition  of  the  Aerotherm  Equilibrium  Surface  Thermo¬ 
chemistry  (EST)  computer  program*.  The  previous  code,  EST  Version 

2,  is  described  in  Reference  1.  The  current  program,  EST  Version 

3,  solves  all  problems  that  earlier  versions  could  solve  and  pro¬ 
vides  additional  computational  capabilities.  The  improved  program 
logic  has  significantly  reduced  machine  run  time  and  decreased  the 
complexity  of  input  preparation. 

The  purpose  of  Volume  I  of  this  user's  manual  is  to 
enable  an  unfamiliar  user  to  effectively  utilize  Version  3  of  the 
Aerotherm  Equilibrium  Surface  Thermochemistry  computer  program. 

To  this  end,  this  document  contains: 

•  a  brief  overview  of  the  purpose  and  capabilities 
of  the  EST  program  (this  section) 

•  a  non-rigorous  discussion  of  the  theoretical 
foundations  of  the  calculation  methods  (Section  2) 

•  a  brief  description  of  the  numerical  solution 
procedures  and  the  FORTRAN  coding  of  these 
procedures  (Section  3) 

•  a  specification  of  the  meanings  and  formats  of 
the  input  data  deck  (Section  4) 

•  some  example  problem  inputs  and  outputs 
(Section  5) . 


This  manual  applies  strictly  to  Version  3  of  the  EST  program  (the 
current  version  as  of  March  1970) .  Users  with  earlier  or  later 
versions  of  the  program  should  contact  the  authors  relative  to 
differences  between  their  program  and  this  manual. 
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Part  II  of  this  manual  contains  the  following  additional 
program  documentations 

•  definitions  for  all  Fortran  variables  used 

•  listings  of  Fortran  IV  source  decks 

•  flow  charts  of  progreim  logic  for  each  Fortran 

routine . 

Three  other  documents  of  potential  value  to  a  reader 
desiring  a  more  detailed  exposition  of  the  theoretical  fundamentals 
of  the  EST  program  are  References  2,  3,  and  4. 

Version  3  of  the  Equilibrium  Surface  Thermochemistry 
program  is  an  extremely  versatile  code  for  calculating  quantities 
of  importance  to  a  broad  variety  of  thermochemical  processes.  The 
thermochemical  processes  treated  may  be  divided  into  two  categories: 
closed  systems  and  open  systems.  Closed  systems  are  those  for 
which  the  relative  amounts  of  each  chemical  element  in  the  system 
is  prespecified.  Open  systems  are  those  for  which  the  relative 
amounts  of  chemical  elements  depend  on  various  mass  transfer  rates 
due,  for  example,  to  boundary  layer  convection  or  solid  surface 
degradation.  The  EST  program  can  treat  both  systems  in  chemical 
equilibrium  and  systems  for  which  sets  of  molecular  species  may  be 
chemically  isolated  from  other  molecules  in  the  system. 

The  variety  of  phenomena  which  may  be  treated  by  the 
EST  program  can  be  illustrated  by  a  hypothetical  example.  Consider 
an  ablative  ma-.erial  (at  the  throat  of  a  rocket  nozzle)  being 
exposed  to  a  solid  propellant  exhaust  stream.  The  situation  is  depic¬ 
ted  in  the  sketch  below.  Encircled  digits  denote  some  of  the  types 
of  calculations  which  may  be  performed  by  the  EST  program.  These 
are  discussed  briefly  below. 

1.  The  thermodynamic  and  chemical  state  at  the 
boundary  layer  edge  may  be  calculated  by 
specifying  the  elemental  composition  (known 
from  the  fuel  and  oxidizer  compositions,  for 
example)  and  two  thermodynamic  state  variables. 
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one  of  which  must  be  pressure  and  the  other 
may  be  static  enthalpy  or  static  temperature. 

This  is  a  closed  system  calculation. 

The  variation  of  enthalpy  with  temperature 
for  the  molecular  composition  of  the  boundary 
layer  edge  gases.  These  calculations  often 
called  the  frozen  edge  table,  are  usually 
required  in  the  evaluation  of  convective 
heat  transfer  to  the  ablative  surface. 

The  thermodynamic  state  of  gas  in  equilibrium 
with  the  surface  of  an  ablating  material  may  be 
calculated  as  a  function  of  pyrolysis  gas  and 
char  mass  rates  normalized  by  the  mass  transfer 
coefficients.  The  effects  of  equal  or  unequal 
diffusion  coefficient  may  be  accurately  accounted 
for  in  the  computation.  These  are  open  system 
calculations  and  the  output  is  usually  in  the 
form  of  tables  for  subsequent  input  to  surface 
energy  balance  calculations. 

The  thermodynamic  state  of  gases  adjacent  to  a 
surface  but  not  in  equilibrium  with  that  surface 
may  be  computed  by  specifying  the  surface  temper¬ 
ature.  Equal  or  unequal  diffusion  may  be  con¬ 
sidered  and  normalized  transpiration  rates  may 
be  allowed.  These  are  also  open  system  calcula¬ 
tions  . 

Equilibrium  Mollier  diagrams  for  any  real  perfect 
gas  (such  as  ablative  material  pyrolysis  gases) 
may  be  computed.  Entropy,  density,  molecular 
weight,  and  either  enthalpy  or  temperature  are 
calculated  as  a  function  of  pressure  and  either 
temperature  or  enthalpy.  Any  atomic,  molecular 
or  ionized  species  including  solids  and  liquids 
Coin  be  considered  in  this  closed  system  calculation. 
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Because  of  the  unusual  versatility  of  the  EST  program, 
the  task  of  instructing  an  unfamiliar  user  in  its  use  is  a  formid¬ 
able  one.  This  manual  endeavors  to  fulfill  this  task  through  a 
brief  summary  of  some  underlying  physical  principles ,  a  specifica¬ 
tion  of  the  input  mechanics,  and  some  sample  problems.  However, 
regardless  of  how  comprehensive  the  manual  and  how  well  checked 
out  the  program,  instances  will  occur  in  which  the  program  fails 
to  solve  a  certain  problem  and  the  reason  for  the  failure  is  not 
apparent  from  the  manual.  These  cases  are  best  resolved  by  direct 
communication  with  the  program  or  manual  authors ,  and  such  communi¬ 
cation  is  encouraged. 
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SECTION  2 
THEORY 


This  section  briefly  summarizes  the  theoretical  founda¬ 
tions  upon  which  the  calculation  methods  of  the  EST  program  are 
based.  An  attempt  has  been  made  to  present  these  fundamentals  in  a 
fashion  most  appropriate  to  a  user  unfamiliar  with  the  code.  To 
this  end,  and  in  the  interest  of  brevity,  explanations  tend  to  be 
nonrigorous  and  perhaps  somewhat  heuristic.  In  discussing  each 
type  of  problem  treated  by  the  EST  program,  the  pertinent  governing 
equations  are  developed  and  it  is  shown  that  these  relations  are 
sufficient  to  determine  the  unknown  quantities  for  that  particular 
problem  type.  No  attempt  is  made  in  this  section  to  discuss  the 
mathematical  methods  employed  to  solve  these  equations.  The  solu¬ 
tion  procedures  are  very  briefly  considered  in  Section  3  which  is 
concerned  primarily  with  the  program  coding. 

The  current  version  of  the  EST  program  is  capable  of 
treating  both  fully  equilibrated  chemical  systems  and  those  for 
which  there  are  mixed  groups  of  isolated  systems  in  individual 
equilibrium.  The  treatment  of  equilibrium  systems  is  highly  gen¬ 
eralized  and  well  checked  out.  The  theoretical  fundamentals 
relative  to  a  variety  of  closed  and  open  equilibrium  systems  are 
discussed  in  Section  2.1.  Techniques  to  obtain  the  composition 
of  systems  containing  isolated  species  which  do  not  react  in  equi¬ 
librium  with  the  main  system  are  described  in  Section  2.2 


* 

Some  judgement  is  required  on  the  part  of  the  program  user  relative 
to  whether  chemical  equilibrium  is  an  appropriate  assumption  or  if 
reaction  kinetics  are  of  importance.  Unfortunately,  there  exists 
no  universally  applicable  criteria  to  aid  the  user  in  making  this 
judgement.  Within  specific  classes  of  problems,  however,  certain 
guidelines  may  be  established  from  experience.  For  example,  for 
the  ablation  of  reinforced  phenolic  materials,  chemical  equilibrium 
is  usually  a  good  assumption  for  surface  temperatures  above  about 
3000°R.  For  graphitic  materials,  reaction  limiting  may  be  important 
up  to  the  onset  of  graphite  sublimation.  A  computer  code  which 
accounts  for  heterogeneous  graphite  reactions  in  solid  propellant 
environments  has  been  developed  under  AF  Contract  F04611-69-C-0081. 
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2 . 1  EQUILIBRIUM  SYSTEMS 

A  variety  of  problem  types ,  subject  to  the  assumption 
of  chemical  equilibrium,  are  discussed  here.  Of  the  many  problem 
options  available  in  the  EST  program,  all  fit  into  one  of  two  cate¬ 
gories:  closed  systems  and  open  systems.  For  closed  systems,  the 
relative  amounts  of  each  chemical  element  is  prespecified.  In 
open  systems,  the  relative  amounts  of  chemical  elements  depend  on 
various  mass  transfer  and  material  degradation  rates.  Closed  sys¬ 
tems  are  discussed  in  Section  2.1.1  and  open  systems  are  discussed 
in  Section  2.1.2. 

2.1.1  Closed  Systems  in  Equilibrium 

The  discussion  of  simple  closed  equilibrium  systems 
in  this  section  is  particularly  significant  because  it  is  here 
that  some  equilibrium  chemistry  relations,  which  are  applicable 
to  all  systems  subsequently  discussed,  are  first  introduced. 

Consider  K  chemical  elements  introduced  into  a  pre¬ 
viously  evacuated  container  in  the  relative  amounts  given  by  N, 

K  • 

In  general,  these  elements  will  interact  to  form  a  number  of 
chemical  species*,  (gas  phase)  and  (condensed  phases). 

If  enough  time  has  elapsed  so  that  thermodynamic  and  chemical 
equilibrium  is  established,  the  thermodynamic  state  of  the  sys¬ 
tem,  including  the  relative  amounts  of  chemical  species  present, 
is  completely  determined  if  two  independent  thermodynamic  var¬ 
iables  are  known**.  Typical  variables  of  this  type  include 
temperature,  pressure,  specific  enthalpy  and  specific  entropy. 


* 

"Chemical  species"  as  used  here  includes  molecular,  atomic, 
ionic,  and  electron  species. 

** 


Duhem's  Theorem. 
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This  condition  may  be  stated  mathematically  by  examining  the  governing 
equations  for  such  a  system,  and  showing  that  the  number  of  independent 
equations  is  equal  to  the  number  of  unknown  quantities . 

Relations  expressing  the  formation  of  the  gaseous  c’  em’ # al 
species  from  the  gaseous  chemical  elements  may  be  written  as  follows: 

K 

Z  \  -  Ni  !i> 

k=l 

Similarly,  formation  of  condensed  phase  species  from  the  gaseous  ele¬ 
ments  is  written: 

K 

Z  nt  <2> 

k=l 


In  the  above,  represents  the  number  of  atoms  of  element  k  in  a 
molecule  of  species  i  (gas)  or  species  Z  (condensed) . 

If  the  gas  phase  species  are  assumed  to  individually 
behave  as  thermally  perfect  gases,  then  the  equilibrium  relation  cor¬ 
responding  to  reaction  (1)  is 

P. 

_ 1 

K 

P 

k=l 


in  pi  -  Z  in  pk  “ in  KPi|T> 

k=l 


Sci 


"  KPi(T) 


(3) 
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where  denotes  partial  pressures  and  K^(T)  is  the  equilibrium  con¬ 
stant  for  the  formation  reaction  (1)  of  species  .  For  each  candidate 
condensed  phase  species 


CkJt  In  Pk  -  in  Kp)!(T) 


where 

=  indicates  the  existance  of  the  condensed  phase  species 
in  equilibrium  with  gas  phase  species,  and 

<  indicates  that  the  condensed  phase  species  will  not 
be  present  in  equilibrium. 

For  each  chemical  element  introduced  into  the  system,  the 
conservation  of  atoms  dictates  that  the  amount  of  any  element  k  in 
the  gas  and  condensed  phases  (regardless  of  molecular  configuration) 
must  sum  to  the  total  amount  of  element  k  in  the  system.  Mathematically, 
this  may  be  written,  for  each  element  k,  as 


Mass  fraction 
of  element  k 
input  to  the 
system 


-  w  E  Vi +  w  E  cmx« 


where  7n  is  a  composite  system  molecular  weight*  defined  by 


xa  (-if •* 9Ss°ofT»r) 

•  .  i  n n 


This  is  the  molecular  weight  appropriate  to  the  ideal  gas  equation 
of  state  if  condensed  phases  are  present. 
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and  where  ^  is  a  mole  fraction  of  condensed  phase  species  l  defined 
as 


x  _  molecules  of  condensed  species  1  (7) 

total  gas  phase  molecules  i 

In  addition,  for  the  gas  phase  species,  there  exists  the  requirement 
that  the  partial  pressures  must  sum  to  the  total  system  pressure 

I 

X  pi  -  f  <» 

i=l 

,  Mixture  thermodynamic  properties,  such  as  specific  en¬ 

thalpy,  are  related  to  the  species  concentrations  by  equations  of 
the  form 


h 


Pihi 


+ 


1 


Vxh 

/Lj  i  l 

«•=! 


(9) 


Consider  now  the  number  of  independent  equations  for  the 
system.  The  number  of  gas  phase  equilibrium  relations  (3)  is  equal 
to  the  number  of  gas  phase  species  I  minus  the  number  of  elements  K 
(because  equations  (3)  are  trivial  when  i=k) .  In  addition,  there  ex¬ 
ists  a  relation  such  as  (4)  for  each  of  the  L  candidate  condensed 
phase  species  in  the  system.  Note  that  the  system  temperature  is 
contained  implicitly  in  equations  (3)  and  (4)  through  the  tempera¬ 
ture  dependence  of  the  equilibrium  constants.  There  are  K  conser¬ 
vation  of  elements  equations  (5) ,  one  for  each  atomic  element  intro¬ 
duced  into  the  system.  The  requirement  that  the  partial  pressures 
sum  to  the  system  pressure  (8)  contributes  one  additional  equation. 
For  any  additional  thermodynamic  properties  of  the  mixture  (enthalpy, 
entropy,  etc.),  there  exists  equations  such  as  (9). 

Consider  next  the  variables  appropriate  to  this  formula¬ 
tion  of  the  problem.  The  relative  concentrations  of  the  I  species 


2-6 


in  the  gas  phase  are  given  by  the  P^s  and  the  amounts  of  the  L 
candidate  condensed  phase  species  are  given  by  X.  (most  or  all  of 

Jv 

which  may  be  zero) .  In  this  formulation,  the  composite  system 
molecular  weight,  %  is  also  a  variable.  There  are  one  each  of 
the  mixture  thermodynamic  variables  T,  P,  h,  s,  etc.  The  number 
of  variables  and  available  independent  equations  may  be  summarized 
as 


VARIABLES 

NO.  OF  SUCH 
VARIABLES 

EQUATION 

NUMBER 

NO.  OF  SUCH 
EQUATIONS 

P. 

1 

1 

(3) 

I  -  K 

L 

(4) 

L 

1 

(5) 

K 

p 

1 

C8) 

1 

T 

1 

h/S/P**** 

n 

of  the  type  (9! 

n 

total 

variables 

I+L+n+3 

total 

equations 

I+L+n+1 

Thus ,  there  are  two  less  equations  than  there  are  variables ,  and  so 
if  two  independent  variables  are  specified  (e.g.,  P  and  T)  in  addi¬ 
tion  to  the  elemental  composition,  then  closure  is  obtained  and  the 
chemical  and  thermodynamic  state  of  the  system  may,  in  principal,  be 
determined. 

The  EST  program  performs  this  determination.  That  is,  to 
determine  the  equilibrium  thermodynamic  and  chemical  state  of  any 
closed  system,  one  needs  only  to  furnish  the  EST  program  with  the 
elemental  composition,  the  candidate  gaseous  and  condensed  phase 
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species,*  and  two  thermodynamic  properties.  One  of  these  properties 
must  be  pressure,  and  the  other  may  be  either  temperature  or  enthal¬ 
py.  Given  this  information,  the  EST  program  will  calculate  and  out¬ 
put  the  mole  fractions  of  each  candidate  species,  the  temperature, 
enthalpy,  entropy,  density,  effective  molecular  weight,,,  equilibrium 
and  frozen  specific  heats,  the  isentropic  exponent,  and  a  few  other 
quantities  of  potential  interest.  Where  appropriate,  these  proper¬ 
ties  are  output  for  the  gas  phase,  condensed  phase,  and  composite 
system.  The  exact  convention  for  inputting  data  to  the  program  is 
specified  in  Section  4 ,  and  some  example  problems  of  this  type  are 
presented  in  Section  5. 

2.1.2  Open  Systems  in  Equilibrium 

The  previous  sections  were  concerned  with  systems  for 
which  the  relative  quantities  of  each  chemical  element  were  pre¬ 
specified.  These  were  termed  closed  systems.  In  many  circumstances 
of  engineering  interest,  however,  the  relative  amounts  of  chemical 
elements  are  not  known  a  prior,; .  Such  systems  are  called  open  sys¬ 
tems.  For  example,  the  gases  adjacent  to  an  ablating  surface  (e.g., 
a  reentry  body  or  a  rocket  nozzle  throat  insert)  constitute  an  open 
system  because  the  elemental  composition  of  these  gases  depend  on 
boundary  layer  transport  events ,  material  pyrolysis  rates ,  material 
char  rates,  etc.  The  EST  program  is  capable  of  treating  this  type 
of  open  system  in  a  particularly  convenient  fashion. 

As  will  be  discussed  in  more  detail  subsequently,  the 
EST  program  relates  the  surface  state  (T  ,  h  ,  etc.)  to  the  ablation 
rate(s) .  The  surface  state  and  ablation  rate(s)  are  also  related  by 
surface  energy  balance  considerations  which  depend  on  the  rate  of 

*The  EST  program  must  be  supplied  certain  basic  thermodynamic  data 
(described  in  Section  4)  for  each  candidate  species  to  be  con¬ 
sidered  in  a  given  system.  This  data  is  contained  in  a  three  card 
set,  one  set  for  each  species.  A  certain  amount  of  judgement  is 
required  on  the  part  of  the  user  relative  to  which  candidate  species 
should  be  included  in  a  given  system.  Frequently  this  judgement 
is  avoided  by  simply  inputting  data  for  all  species  containing  com¬ 
binations  of  the  input  elements. 
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heat  conduction  into  the  solid.  Thus,  the  particular  ablation  rate, 
surface  temperature,  etc.,  is  calculated  by  combining  these  two 
relations:  surface  thermochemistry  +  in-depth  heat  conduction.  In 
addition,  the  boundary  layer  heat  and  mass  transfer  coefficients 
must,  of  course,  be  specified.  Computationally,  the  surface  thermo¬ 
chemistry  relations  are  usually  generated  by  the  EST  program  in  the 
form  of  tables  (card  sets)  for  subsequent  input  to  transient  heat 
conduction  and  ablation  programs  such  as  the  CMA  (Charring  Material 
Ablation  -  Reference  5)  program,  the  SCRIMP  (Silica  Carbon  Reactions 
Including  Melting  Phenomena  -  Reference  6)  program,  or  the  ASTHMA 
(Axisymmetric  Transient  Heating  and  Material  Ablation  -  Reference  7) 
program.  The  information  flow  for  a  complete  transient  ablation 
prediction  may  be  illustrated  diagramatically  as: 


surface  temperature  and  recession 
history,  in-depth  temperature 
distribution  history,  etc. 


The  preparation  of  the  input  for  these  ablation  calcula¬ 
tions  is  discussed  in  Section  4  for  the  EST  program  and  in  References 
5,  6,  and  7  for  CMA,  SCRIMP,  and  ASTHMA  programs,  respectively.  The 
fundamentals  underlying  open  system  thermochemistry  are  briefly 
summarized  here.  Open  system  calculations,  subject  to  a  few  simpli¬ 
fying  assumptions  to  more  clearly  illustrate  the  basic  theory,  are 
discussed  first  in  Section  2. 1.2.1.  These  considerations  are  exten¬ 
ded  to  more  general  systems  in  Section  2. 1.2. 2. 
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2. 1.2. I  Open  Systems  in  Equilibrium  -  Simplified  Case 

The  basic  theory  underlying  the  treatment  of  open  systems 
may  best  be  illustrated  by  examining  the  equations  expressing  the 
conservation  of  chemical  elements  and  energy  at  the  ablating  surface. 

If  the  boundary  layer  is  characterized  by  equal  diffusion  coefficients, 
unrty  Prandtl  number,  and  unity  Lewis  number,  and  if  no  material  is 
removed  from  the  surface  in  a  condensed  phase  (i.e.,  no  mechanical 
erosion  or  liquid  layer  removal) ,  then  these  equations  take  on  a 
particularly  simple  form  for  equilibrium  systems.  This  simplified 
situation  will  be  considered  in  this  section  and  these  considerations 
will  be  extended  to  more  generalized  cases  in  Section  2. 1.2. 2. 

Consider  first  the  fluxes  of  chemical  elements  (k) 
entering  and  leading  a  control  surface  affixed  to  the  ablating 
surface.  The  solid  material  may  be  visualized  as  moving  into  this 
surface  at  a  rate  s.  If  it  is  assumed  that  no  material  is  being 
removed  in  a  condensed  phase,  then  the  surface  and  the  fluxes  of 
the  k  chemical  element  may  be  illustrated  as  follows: 


boundary 
layer  edge 


ablating 

surface 


s 
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Terms  superscripted  by  a  tilde  (~)  represent  the  total  mass  fraction 
or  flux  of  element  k,  independent  of  molecular  configuration.  Thus 


vKi 

i  w 


£ 

■i 


'k.  Ji 
1  w 


(10) 


(11) 


where  k  pertains  to  element  k,  i  pertains  to  species  i,  and  is 

the  mass  fraction  of  element  k  in  species  i.  Fluxes  of  element  k  away 

from  the  surface  consist  of  boundary  layer  diffusion  and  gross  motion 
of  the  fluid  adjacent  to  the  surface  due  to  injection  fluxes  m  (pyrol- 
ysis  gas  rate)  and  mc  (char  rate) . 

From  the  above  sketch,  it  is  seen  that  conservation  of 
chemical  elements  requires  that 

3k  +  (pvl»“k  -  Vk  +  “A  1121 

w  w  ^  g  c 

Summing  Equation  (12)  over  all  elements  k  yields  the  total  mass  con¬ 
tinuity  equation  (for  the  case  when  there  is  no  condensed  phase  mate- 

rj  al  removal) 

(pV)w  =  ^g  +  “c  (13) 

An  important  fundamental  of  the  present  mathematical  modeling  of  the 
ablation  process  is  the  expression  of  the  diffusive  heat  and  mass 
fluxes  in  terms  of  a  transfer  coefficient  formulation.  This  formula¬ 
tion  will  be  discussed  briefly  in  the  following  paragraphs,  and  more 
detailed  treatments  are  given  in  References  8,  9,  and  10. 

Heat  and  Mass  Transfer  Coefficients  -  Simplified  Case 

For  low  speed  flow  of  an  incompressible,  constant 

property,  non-reacting  fluid,  the  boundary  layer 
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energy  equation  may  be  written  in  the  form 


u 


ot 


ill 

3y2 


(14) 


and  solutions  to  this  equation  have  historically  been 
correlated  in  terms  of  a  heat  transfer  coefficient  h 
where 


q  =  h(Te  -  Tw)  (15) 

For  high  speed  chemically  reacting  boundary  layers  (as 
are  of  interest  in  ablation  problems) ,  the  energy  equa¬ 
tion  can  be  written  in  the  following  form  if  diffusion 
coefficients  are  equal  and  if  the  Prandtl  and  Lewis 
numbers  are  unity  (e.g.,  Reference  8) 

3H  x  3H  _  3  /  3H\ 

pu  3S  +  pv  17  3y  V  ^7  (  6) 

where  H  is  the  total  (sensible  +  chemical  +  u2/2)  enthalpy. 
By  analogy  to  Equations  (14)  and  (15) ,  solutions  to  (16) 
are  conveniently  expressed  in  terms  of  a  dimensionless 
heat  transfer  coefficient  C  where 

fi 

q  -  »e%CH(He  -  Hw>  (17) 

Corresponding  formulations  for  non-unity  Prandtl  and 
Lewis  number  cases  will  be  considered  in  Section  2. 1.2. 2. 

Consider  now  the  equation  for  mass  transfer  diffusion  in 
the  boundary  layer.  If  the  diffusion  coefficients  for 
all  chemical  species  are  equal  (Pick's  Law) ,  then  an 
appropriate  summation  (Reference  8)  of  the  boundary 
layer  species  conservation  equation  yields  the 
following  equation  for  each  element  k 
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l 

3x  +  pv  3y  _  3y 


(18) 


Again,  by  analogy,  solutions  to  this  equation  may  be  cor¬ 
related  in  terms  of  a  dimensionless  mass  transfer  coeffi¬ 
cient  CM  where 

-Sk  "  peueCM^Kk  -  h  >  <ls> 
w  e  w 


And  for  the  unity  Prandtl  and  Lewis  numbers  case  under 
consideration,  the  similarity  of  Equations  (16)  and  (18) 
indicates  that  if  the  corresponding  boundary  conditions 
are  also  similar,  then  CM  =  C„. 


Utilizing  this  transfer  coefficient  formulation  (19)  for 
the  diffusion  flux  in  the  elemental  balance  Equation  (12)  yields 


•>eueCM<\  -  \e>  +  (pv>„\  -  Vk  + 


(20) 


defining  dimensionless  pyrolysis  and  char  rates  as 


m 

B'  =  - 2— _ 

^  peueCM 


m 


B '  =  - 

c  p  u  C. 


and  B '  s 


(pv) 


w 


e  e  M 


peuecM 


(21) 


and  solving  (20)  for  the  total  mass  fraction  of  element  k  at  the  wall 
yields  (for  equal  diffusion  coefficients,  Pr  =  Le  =  l.o,  and  no  con¬ 
densed  phase  material  removal) 


Vk  +  B;Kk 

g _  c 

1  +  B' 


+  K,. 


(22) 


Given  the  relative  amounts  of  chemical  elements  specified 
by  (22)  ,  the  chemical  and  thermodynamic  state  of  the  gases  adjacent 
to  the  ablating  surface  may  be  calculated  from  equilibrium  relations 
similar  to  those  discussed  in  Section  2.1.1. 


r 


i 

t 


}■ 

t 
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Since  the  gases  are  in  equilibrium  with  the  ablating  surface 

K 

-  £  Ckt  ln  h  *  ln  V(T) 

k=l 

if  £  represents  the  surface  species,  and 

K 

-  £  h ln  pk  ‘  v(T> 

k=l 


(23) 


(24) 


for  all  other  candidate  condensed  phase  species.  In  the  present  for¬ 
mulation,  (24)  may  be  thought  of  as  being  used  to  identify  surface 
species  £  for  which  (23)  applies.  The  equilibrium  relations  for  gas 
phase  species  is  again  (see  Section  2.1.1)  of  the  form 

K 

tn  pi  -  Z  In  pk  -  in  KPi (t>  ,25> 

k=l 


and  the  partial  pressures  must  obey  the  relation 


I 

£  Pi  ■  r  (26) 

i=l 

The  elemental  mass  fractions  adjacent  to  the  surface,  Kj^,  of  (22) 
are  related  to  the  species  partial  pressures,  P^,  by  relations  such 
as  (5)  (except  that  no  condensed  phase  is  being  considered  here) 


(27) 


Thus,  if  P  is  known  and  and  are  specified  (these  may  be  varied 
parametrically  as  will  be  discussed  subsequently)  ,  and  if  Tw  and 
are  unknowns,  then  the  number  of  unknowns  and  equations  available 
may  be  summarized  as 


2-14 


UNKNOWNS 

NO.  OF  SUCH 
UNKNOWNS 

EQUATION 

NO. 

NO.  OF  SUCH 
EQUATIONS  AVAILABLE 

P. 

1 

I 

(25) 

I-K 

K 

(23) 

1 

1 

(22) 

K 

T 

w 

1 

(26) 

1 

(27) 

K 

Total 

Unknowns 

I+K+2 

Total 

Equations 

I+K+2 

Thus,  closure  is  obtained  and,  in  principal,  the  tempera¬ 
ture  and  chemical  composition  of  the  gases  adjacent  to  the  surface 

may  be  determined  if  P,  B',  and  B'  are  specified.  From  the  pressure, 

g  c 

temperature,  and  chemical  composition,  the  calculation  of  other 
thermodynamic  properties  (enthalpy,  etc.)  is  straightforward. 


The  EST  program  has  the  capability  to  determine  the  chem¬ 
ical  and  thermodynamic  state  of  the  gases  adjacent  to  an  ablating  sur¬ 
face  in  a  fashion  similar  to  that  discussed  here.  Indeed,  the  EST 
program  is  capable  of  treating  more  complex  ablation  situations, 
and  these  are  discussed  in  Section  2. 1.2. 2.  In  general,  the 

pressure  P  and  a  matrix  of  B'  and  B'  values  are  input.  This  is 

g  c 

usually  most  convenient  since  Tw  is  usually  single  valued  with  re¬ 
spect  to  B^,  and  Bg  as  illustrated  by  a  typical  example 


B'  and  P  constant 

g 


Clearly,  calculating  from  input  Tw  would  not  be  appropriate  for 
this  case .  For  other  cases  such  as 
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the  program  possesses  some  special  "fill-in"  options  (discussed  in 
Section  4) .  For  open  systems ,  the  EST  output  consists  of  card  sets 
(in  addition  to  regular  printed  output)  referred  to  as  surface 
thermochemistry  tables .  Each  card  contains  information  relative  to 
one  surface  state  (i.e.,  one  P,  T ,  B',  B',...,  combination).  As 

w  y  C 

previously  discussed,  this  surface  thermochemistry  data  alone  is 
insufficient  to  solve  an  ablation  problem;  it  must  be  combined  with 
in-depth  heat  conduction  and  surface  energy  balance  considerations. 

In  order  to  more  clearly  demonstrate  this  interdependence,  the 
relation  of  the  surface  thermochemistry  to  the  surface  energy  balance 
is  discussed  briefly  in  the  following  paragraphs. 


Relation  to  Surface  Energy  Balance  -  Consider  the  fluxes 

of  energy  entering  and  leaving  a  control  surface  affixed 

to  the  ablating  surface.  For  the  unity  Prandtl  number, 

unity  Lewis  number,  no  condensed  phase  removal  ablation 

case  being  considered,  these  fluxes  may  be  illustrated  as 

boundary 
layer  edge 

q  ( pv) H  q  ,  .  aeT  4 

^conv  w  ^rad  in  w 


ablating 

surface 


'cond 


m  H 

g  g 


m  H 
c  c 
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where: 

qconv  is  the  ener9Y  convected  to  the  surface  as  a 

result  of  boundary  layer  transport  events. 

This  term  includes  effects  of  both  heat  con¬ 
duction  due  to  a  temperature  gradient  in  the 
gas  adjacent  to  the  surface,  and  energy  trans¬ 
port  due  to  endothermic  or  exothermic  chemi¬ 
cal  reactions  at  the  surface.  For  the  equal 
diffusion  coefficient,  Pr  =  Le  <=  1.0  case 
under  consideration  here,  the  previously  dis¬ 
cussed  transfer  coefficient  form  of  this  term 
is  (17) 

^conv  "  "eW’e  "  V  •  <l7> 

(pv)wHw  is  the  energy  flux  leaving  the  surface  with 
the  "blowing"  flux  (pv)^, 

q^_.  •!„  is  the  incident  radiation  flux, 

rad  m 

aeT^1*  is  the  surface  re-radiation  flux, 

q  ,  is  the  heat  flux  conducted  into  the  solid 

cond  9T\ 

material,  q,  d  =  "  *  ^ 

J/w, solid 

m  H  and 
9  9 

m  H  are  the  energy  fluxes  entering  the  surface 

c  c 

associated  with  the  pyrolysis  gas  and  char 
respectively. 

Conservation  of  energy  at  the  surface  requires  that 

peueCH(He  "  V  +  “gHg  +  ®oHc  ”  (pv)wHw  +  ^rad  in  “ 

oeV  =  ^cond 


or,  utilizing  (21) 
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For  most  problems,  He  is  known,  PeueCH  is  independently 

determined  from  boundary  layer  calculations,  e  is  known, 

qrad  is  known  (frequently  zero) ,  and  Hg  and  Hc  are  known 

functions  of  T  .  Additionally,  B'  and  H.,  are  related  to 
w  c  w 

Tw  through  the  surface  thermochemistry  tables  generated 
by  the  EST  program  as  previously  discussed.  Based  on 
considerations  thus  far,  all  quantities  in  (28)  except 
qcond  and  Bg  are  determined  for  a  given  The  in-depth 

heat  conduction  term,  qcond>  depends  only  on  previous  heat 
conduction  events  and  the  current  Tw.  The  pyrolysis  gas 
rate,  PV  ,  is  given  by  a  decomposition  law  (Reference  9) 
which  depends  on  in-depth  temperature  which,  through  heat 
conduction  considerations,  depends  on  Tw.  Thus,  it  is 
seen  that  in  general  there  exists  only  one  value  of  Tw 
which  satisfies  (28) .  This  Tw  is  usually  determined  in 
an  iterative  fashion  by  a  heat  conduction  and  ablation 
program  (e.g..  References  6  or  7)  which  utilizes  the  sur¬ 
face  thermochemistry  tables  as  boundary  conditions.  Once 
Tw  is  known  for  a  given  time  in  a  transient  solution,  B^, 
B',  surface  recession  rate,  in-depth  temperature  distri- 
butions,  and  other  quantities  of  interest  follow  directly. 

The  above  has  been  a  brief  discussion  indicating  the  in¬ 
terdependence  of  the  surface  thermochemistry  solution  and  the  in-depth 
heat  conduction  solution  for  a  somewhat  simplified  case.  It  should 
be  pointed  out  that  similar  considerations  apply  for  situations  when 
unequal  species  diffusion  effects  are  significant  and  when  Pr  ^  0  ^ 
1.0.  The  EST  program  is  able  to  treat  these  effects  and  details 
relative  to  the  treatment  of  them  will  be  considered  in  Section 
2. 1.2. 2 
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2.1.2. 2  Open  Systems  in  Equilibrium  -  Nonunity  Prandtl  and 

Lewis  Numbers ,  Unequal  Diffusion  Coefficients 

The  discussion  of  the  previous  section  was  limited 
to  open  systems  with  unity  Prandtl  and  Lewis  numbers,  equal 
species  diffusion  coefficients,  and  no  removal  of  surface  mat¬ 
erial  in  the  condensed  phase.  These  simplifications  were  made 
in  order  to  render  the  basic  theory  easier  to  explain.  While 
this  simple  model  is  reasonably  accurate  for  many  ablation  sit¬ 
uations,  these  assumptions  are  inappropriate  for  others.  The 
equal  diffusion  coefficient,  Pr  =  Le  *  1.0  simplifications  per¬ 
tain  to  boundary  layer  mass,  energy,  and  momentum  transport 
events,  and  the  effects  of  the  relaxation  of  these  assumptions 
on  the  problem  formulation  and  solution  procedure  will  be  briefly 
discussed  in  this  section. 

Pr  ^  Le  ^  1.0  -  For  nonunity  Prandtl  and  Lewis  num¬ 
bers,  the  transfer  coefficient  formulation  for  the  boundary 
layer  energy  flux  is  not  as  straightforward  as  that  discussed  in 
the  previous  section.  This  is  because  the  boundary  layer  energy 
equation  is  no  longer  of  the  similar  form  and  thus  a  transfer 
coefficient  formulation  cannot  be  justified  purely  by  analogy. 

A  detailed  discussion  of  boundary  layer  transport  models  for  the 
Pr  f  Le  1  1.0  case  is  beyond  the  scope  of  this  manual  (  see,  e.g., 
References  9  and  10)  and  only  a  few  of  the  results  as  they  relate 
to  film  transfer  coefficient  formulations  will  be  discussed  here. 

When  the  Prandtl  number  is  not  unity,  the  viscous 
dissipation  and  heat  conduction  terms  in  the  boundary  layer  energy 
equation  cannot  be  combined  thus  rendering  the  equation  inhomo¬ 
geneous.  Solutions  to  this  equation  indicate  that  the  driving 
potential  in  the  transfer  coefficient  expression  for  the  surface 
heat  flux  should  be  defined  in  terms  of  a  recovery  enthalpy  (e.g., 
Reference  11)  in  place  of  the  actual  boundary  layer  edge  enthalpy. 
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Thus,  nonunity  Prandtl  number  has  no  influence  on  EST  calculations 
and  the  primary  practical  implication  is  that  the  recovery  enthalpy, 
rather  than  the  boundary  edge  enthalpy,  should  be  input  to  the  heat 
conduction  and  ablation  solution. 


When  the  Lewis  number  is  not  unity,  the  terms  in  the 
boundary  layer  energy  equation  representing  energy  transfer  by  heat 
conduction  and  chemical  species  diffusion  cannot  be  combined,  again 
rendering  the  equation  inhomogeneous.  The  energy  flux  to  the  surface 
is  given  by 


hj°  pD 


(29) 


where  the  first  term  characterizes  the  heat  conducted  to  the  surface 
as  a  result  of  the  temperature  gradient  in  the  gas  adjacent  to  the 
surface,  and  the  second  term  represents  the  effect  of  endothermic 
and  exothermic  chemical  reactions  at  the  surface.  The  appropriate 
transfer  coefficient  form  of  (29)  is  not  firmly  established  at  this 
time.  However,  Reference  9  suggests  the  form 
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(30) 


and  this  is  the  form  employed  in  Aerotherm  ablation  programs.  In 
(30) ,  the  driving  potential  in  the  first  term  is  the  recovery  en¬ 
thalpy  at  the  boundary  layer  edge  minus  the  enthalpy  of  the  boundary 
layer  edge  gases  frozen  at  the  edge  composition  and  at  the  surface 
temperature,  and  h^w  represent  the  enthalpy  of  chemical  species  i  with 
respect  to  a  base  temperature  equal  to  the  surface  temperature.  It  can 
be  shown  (Reference  9)  that  for  Le  =  1.0  and  CM  *  CH,  (30)  collapses 
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to  (17)  as  expected.  However,  where  Le  ^  1.0,  the  heat  and  mass 
transfer  coefficients  are  generally  unequal  and  a  correlation  fre¬ 
quently  employed  (Reference  12)  is 


C 


M 


CH 


(31) 


It  is  apparent  from  (30)  that,  in  addition  to  surface  thermochemical 
data,  boundary  layer  edge  thermochemical  data  are  also  required  when 
Le  f  1.0  in  order  to  specify  the  quantity 

^  T 

hw  frozen  =  L,  Hehi 
edge  gas  i=l 

This  quantity  is  also  calculated  by  the  EST  program  for  subsequent 
input  to  heat  conduction  and  ablation  programs .  The  calculation  of 
these  frozen  edge  tables  is  usually  accomplished  by  performing  a 
closed  system  equilibrium  calculation  at  the  boundary  layer  edge  temp¬ 
erature  and  elemental  composition  (as  discussed  in  Section  2.1.1)  and 
then  performing  closed  system  calculations  frozen  at  this  composition 
for  an  array  of  temperatures  spanning  the  expected  wall  temperature 
range.  The  program  input  for  frozen  composition  calculations  is  dis¬ 
cussed  in  Section  4.  Thus,  for  problems  with  CM  ^  C„,  frozen  edge 
tables,  as  well  as  surface  thermochemistry  tables,  are  prepared  by 
the  EST  program  for  input  to  the  CMA  (Reference  5) ,  ASTHMA  (Reference 
7) ,  or  related  heat  conduction  and  ablation  programs. 

Unequal  Diffusion  Coefficients  -  A  significant  simplifi¬ 
cation  of  the  boundary  layer  energy  and  mass  diffusion  equations  re¬ 
sults  if  all  binary  diffusion  coefficients  for  a  given  species  i  are 
assumed  to  be  equal  (e.g.,  Reference  11)  and  all  considerations  up 
to  this  point  have  been  predicated  on  this  assumption.  However,  for 
many  chemical  systems  of  interest  in  ablation  problems  (e.g.,  when 
there  is  a  significant  difference  between  the  molecular  weights  of 
the  major  species  present)  this  assumption  is  a  severe  compromise 
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with  reality.  Based  upon  an  accurate  approximation  for  binary  dif¬ 
fusion  coefficients,  Reference  3  presents  simplified  equations  for 
a  multicomponent  boundary  layer  with  unequal  diffusion  coefficients 
for  all  species.  The  application  of  similarity  arguments  to  these 
equations  suggests  the  following  form  for  the  transfer  coefficient 
formulation  for  the  diffusion  flux  of  element  k  at  the  surface 

^  -  wM(*k  -  h  >  <33) 

w  we 


In  (33),  Z*  is,  in  effect,  a  weighted  average  of  the  mole  and  mass 
fractions  of  element  k.  The  Z£  are  defined  by 
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(y  ~  2/3,  see 
Reference  2) 
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where  the  F^  are  diffusion  factors  defined  by  the  following  relation 
for  the  binary  diffusion  coefficients 


D 

F.  FT 
i  3 


(37) 


where  D  is  a  constant  for  a  given  temperature  and  pressure  and  the 
Fi  depend  weakly  on  temperature.  The  ^ j  must  obey  (37)  in  order 
for  the  boundary  layer  species  diffusion  equations  to  reduce  to  a 
form  from  which  (33)  can  be  inferred  by  similarity  arguments.  Ref¬ 
erence  3  demonstrates  that  the  binary  diffusion  coefficients  for  a 
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variety  of  chemical  systems  are  accurately  correlated  by  (37).  This 
reference  also  shows  thu  reasonably  good  correlation  equation  for 


the  is 


Fi  * 


^ref 


where 


^ref  ^  ^3.4  and  e  «  0.431  (38) 


when  D  is  taken  as  the  self-diffusion  coefficient  of  0 2* 

Consideration  of  unequal  diffusion  coefficients  affects  the 
surface  elemental  balance  relationships  which  are  needed  to  determine 
the  equilibrium  thermochemical  state  at  the  surface.  Substituting 
(33)  into  (12)  yields 


peueCM(i£  "  KJ  +  <l>viwKk  -  mgKk  +  moKk 
we  w  3  g  c 


and  the  unknowns  here  are  Kk  and  Zk  ,  each  of  which  may  be  expressed 

W  v7 

in  terms  of  the  species  partial  pressures 
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where  is  the  mean  molecular  weight  of  the  gas  phase  and  F  is  a 
mean  F?  defined  as 


F  = 


(42) 


Substituting  these  expressions  into  (39)  and  utilizing  (21)  yields 
an  expression  for  the  species  partial  pressures  at  the  surface  in 
terms  of  quantities  at  the  boundary  layer  edge  and  in  the  material 


(43) 


Note  that  (43)  reduces  to  (22)  when  the  diffusion  coefficients  are 
equal . 


When  performing  unequal  diffusion  coefficient  open  system 
calculations,  the  EST  program  utilizes  (43)  rather  than  (22)  as  the 
elemental  mass  balance  equations.  Other  than  this,  the  solution  phi¬ 
losophy  is  essentially  as  discussed  in  Section  2. 1.2.1.  The  diffusion 
factors  utilized  in  the  EST  program  may  be  calculated  in  three  ways, 
at  the  user's  option 

a.  diffusion  factors  F^  may  be  input  individually  for 
each  species  i 

b.  diffusion  factors  may  be  calculated  according  to 
(38)  with  the  user  specifying  the  reference  molecu¬ 
lar  weight,  7!\re and  the  exponent  e 

c.  if  neither  of  the  above  options  are  utilized 
the  program  will  calculate  F^  according  to  (38) 
with  7i\  f  =  23.4  and  e  =  0.431. 
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The  actual  program  input  for  these  alternatives  will  be  discussed 
in  Section  4. 

For  unequal  diffusion  coefficients,  the  transfer  co¬ 
efficient  formulation  for  the  surface  energy  flux  has  the  follow¬ 
ing  form  (Reference  9) 
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Note  that  for  equal  diffusion  coefficients,  Z*  =  K.  and  (44)  reduces 
to  (30)  as  expected.  Consistent  with  (44),  for  unequal  diffusion 
coefficient  problems,  the  surface  thermochemistry  tables  prepared 
by  the  EST  program  contain  the  quantity 


£ 


in  addition  to  the  quantities  previously  discussed.  Note  again,  that 
for  equal  diffusion  coefficients 
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Similarly,  the  previously  discussed  frozen  edge  table  must  contain 
the  quantity 
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in  addition  to 


h  frozen 
w  edge  gas 

The  EST  program  input  arrangement  for  the  preparation  of  these 
surface  and  edge  thermochemistry  tables ,  for  equal  or  unequal  dif¬ 
fusion  coefficients,  is  specified  in  Section  4. 

2.2  NON-FULLY  EQUILIBRATED  SYSTEMS 

To  calculate  the  equilibrium  state  of  a  chemical  system, 
information  relative  to  all  chemical  reactions  is  not  needed.  This 
fact  permits  a  significant  simplification  in  the  problem  formulation, 
and  those  options  of  the  EST  program  discussed  thus  far  take  advan¬ 
tage  of  these  simplifications .  For  some  systems  of  interest,  how¬ 
ever,  the  effects  of  reaction  kinetics  may  not  be  neglected.  A 
general  solution  of  complex  problems  for  which  reaction  kinetics 
effects  are  important  is  potentially  difficult  for  at  least  two  rea¬ 
sons:  a)  there  are  significant  computational  and  bookkeeping  problems 
associated  with  the  analytical  treatment  of  mixed  equilibrium  and  non¬ 
equilibrium  systems,  and  b)  for  many  systems  of  engineering  interest, 
the  rate  controlling  reactions  are  not  well  known  and/or  rate  con¬ 
stants  for  these  reactions  are  unavailable.  Difficulties  of  these 
types  frequently  preempt  any  exact  analytical  treatment.  In  these 
cases ,  various  approximate  treatments  can  often  yield  useful  informa¬ 
tion.  These  approximate  treatments  essentially  consist  of  not  allow¬ 
ing  certain  species ,  or  groups  of  species ,  to  react  with  an  otherwise 
equilibrated  system.  A  few  examples  of  this  sort  of  approximate 
treatment  will  be  discussed  in  the  following  paragraphs. 

Completely  Frozen  Systems  -  The  EST  program  contains  a 
provision  for  performing  thermodynamic  state  calculations  for  sys¬ 
tems  with  frozen  chemical  compositions.  The  calculations  are  carried 
out  by  performing  an  equilibrium  calculation  at  one  temperature  and 
pressure,  freezing  the  chemical  composition,  and  then  performing 
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thermodynamic  state  calculations  at  other  temperatures  and  pres¬ 
sures  but  at  this  frozen  chemical  composition.  This  type  of  calcu¬ 
lation  was  previously  discussed  in  reference  to  frozen  boundary 
layer  edge  tables  (in  Section  2. 1.2. 2).  Program  input  details 
relative  to  this  type  of  calculation  are  discussed  in  Section  4. 

Removal  of  Selected  Species  from  the  System  -  the 
simplest,  and  perhaps  the  crudest,  way  to  approximate  certain  kin¬ 
etic  effects  is  to  simply  remove  from  the  equilibrium  system  those 
species  whose  formation  is,  in  reality,  suppressed  by  reaction  kin¬ 
etic  effects.  Computationally,  these  species  are  removed  merely 
by  removing  the  thermochemical  data  cards  for  that  molecule  from 
the  species  thermochemical  data  deck.  See  Section  4. 

Isolated  Species  -  Certain  chemical  species  may  be  fro¬ 
zen  individually  at  any  given  concentration.  This  treatment  is 
sometimes  useful  for  species  which  are  relatively  nonreactive  be¬ 
cause  of  chemical  kinetic  effects.  For  example,  consider  the  flow  of 
of  a  C-O-H  gas  over  a  relatively  low  temperature  ablating  carbon 
surface.  Isolation  of  H20  at  its  boundary  layer  edge  concentration 
would  provide  a  more  realistic  surface  thermochemistry  table  than  if 
full  equilibrium  were  assumed  (since  the  water-gas  reaction  at  the 
surface  is  relatively  slow  at  lower  temperatures) .  In  the  EST  pro¬ 
gram,  this  isolation  of  species  is  achieved  by  treating  that  parti¬ 
cular  species  as  an  atomic  element,  an  element  with  a  fictitious 
atomic  number  so  that  it  cannot  react  chemically  with  any  other  mem¬ 
bers  of  the  chemical  systems.  Computationally,  this  requires  speci¬ 
fication  of  the  input  of  the  appropriate  quantity  of  this  "element" 
(accompanied  by  appropriate  reductions  in  the  amounts  of  other  actual 
elements) ,  and  the  alteration  of  the  species  thermochemistry  data 
card  set  for  that  species  so  that  it  contains  only  one  "atom"  of  the 
fictitious  "element". 

Isolated  Subgroups  of  Elements  -  In  some  practical  appli¬ 
cations  ,  there  exists  a  subgroup  of  elements  which  are  in  equilibrium 
with  each  other,  but  the  subgroup  is  not  in  chemical  equilibrium  with 
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the  system  as  a  whole.  For  example,  there  are  indications  (Reference 
13)  that  for  some  ablative  plastics,  the  pyrolysis  gases  do  not  come 
into  chemical  equilibrium  with  the  boundary  layer  gases.  This  effect 
may  be  studied  by  performing  a  solution  with  the  pyrolysis  gases  iso¬ 
lated  from  the  rest  of  the  system.  The  isolation  of  a  subgroup  of 
elements  can  be  achieved  by  considering  these  elements  to  be  differ¬ 
ent  "elements"  with  properties  the  same  as,  but  with  atomic  numbers 
different  from,  the  corresponding  actual  elements  (e.g.,  add  100  to 
the  atomic  numbers  of  the  elements  in  the  isolated  subgroup) .  For 
this  type  of  calculation,  the  subgroup  to  be  isolated  is  considered 
to  be  one  component,  and  the  composition  of  this  component  is  speci¬ 
fied  as  the  actual  relative  elemental  composition  of  the  subgroup  to 
be  isolated,  but  in  terms  of  the  altered  elements.  In  order  that 
these  altered  elements  may  equilibrate  with  each  other,  the  species 
thermochemical  data  card  set  must  contain  a  set  of  species  which  may 
be  formed  from  these  altered  elements,  i.e.,  species  with  the  same 
atomic  number  convention. 

Low  Temperature  Surface  Equilibrium  Suppression  -  As 
previously  discussed,  equilibrium  is  an  increasingly  poor  assumption 
at  low  surface  temperatures  (e.g.,  below  2000°R)  in  some  ablation 
problems.  However,  when  generating  surface  thermochemistry  tables 
with  the  EST  program  (Section  2. 1.2. 2),  lower  temperature  ranges 
are  usually  only  of  minor  interest  (e.g.,  during  a  very  short  ini¬ 
tial  period  of  a  transient  ablation  problem) .  Nontheless ,  the  sur¬ 
face  tables  must  frequently  extend  to  these  lower  temperatures  in 
order  to  be  compatible  with  the  subsequent  heat  conduction  and  abla¬ 
tion  solution.  For  these  cases,  the  EST  program  contains  an  option 
by  which,  for  each  pyrolysis  gas  rate  and  at  the  lowest  specified 
char  rate,  surface  equilibrium  will  be  suppressed  and  a  series  of 
specified  temperature  solutions  will  be  performed  (at  500 °K  intervals) 
to  fill  the  table  out  down  to  500 °K.  For  many  cases,  surface  thermo¬ 
chemistry  tables  generated  in  this  fashion  are  closer  to  reality  than 
full  equilibrium  solutions  down  to  500 °K.  The  input  convention  for 
implementing  this  option  is  discussed  in  Section  4. 
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SECTION  3 
CODING 

This  section  presents  some  details  relative  to  the 
FORTRAN  coding  of  the  EST  program.  A  brief  description  of  the 
numerical  solution  procedure,  which  includes  a  discussion  of  the 
"base  species"  concept,  is  presented  first  in  Section  3.1.  Each 
subroutine  making  up  the  EST  program  is  identified  and  briefly 
discussed  in  Section  3.2.  Definitions  of  several  diffusion  factor 
dependent  properties  appearing  on  the  program  output  are  given  in 

Section  3.3.  Some  miscellaneous  program  details  are  documented  in 
Section  3.4.  These  include  storage  requirements,  tape  requirements, 
operator  controls,  dumps  and  diagnostics,  and  run  time. 

3 . 1  COMPUTATIONAL  PROCEDURE 

The  computational  procedure  employed  in  the  EST  program 
to  solve  the  equations  set  forth  in  Section  2  is  briefly  discussed 
here.  Considerably  greater  detail  is  presented  in  Reference  4  and 
in  particular  Table  1  of  that -reference. 

3.1.1  Basic  Solution  Technique 

The  basic  solution  technique  may  be  illustrated  by  con¬ 
sidering,  for  example,  an  open  system  with  unequal  diffusion 
coefficients  (i.e.,  as  discussed  in  Section  2. 1.2. 2).  For  this 
system,  the  basic  equations  defining  the  problem  are  the  elemental 
conservation  equations  (43)  ,  the  total  pressure  equation  (26)  ,  the 
reaction  equilibrium  equations  (25)  ,  and  one  heterogeneous  vapor 
pressure  relation  (23).  The  table  of  Section  2. 1.2.1  shows  that 
there  are  as  many  knowns  as  unknowns  in  these  equations  so  closure 
is  obtained. 


The  number  of  unknowns  could  immediately  be  reduced  by 
I-K  through  the  direct  substitution  of  (48),  as  solved  for  P^,  into 
the  other  relations „  It  proves  advantageous,  however,  to  continue 
to  treat  the  full  set  of  equations,  and  to  subsequently  utilize  this 
substitution  during  the  iterative  convergence  procedure.  The  solu¬ 
tion  of  these  simultaneous  nonlinear  algebraic  equations  is  based  on 
Newton-Raphson  iteration.  Since  this  procedure  is  accelerated  by 
casting  the  equations  into  a  more  linear  form  (via  transformations, 
substitutions,  etc.)  it  is  well  to  examine  the  set  of  equations 
above.  With  the  boundary  layer  edge,  char  and  pyrolysis  gas  compo¬ 
sition  given  as  well  as  the  B' ,  (46)  and  (47)  are  linear  relations 
between  the  P^  and  ^  providing  that  F  is  reasonably  constant.  In 
contrast,  (48)  and  (49)  are  li.^ar  relations  between  the  5-n  P^, 

£n  P^  and  S.n  Kp^ ,  the  latter  variable  being  approximately  linear  in 
1/T. 

The  EST  program  takes  advantage  of  this  situation  by 
treating  those  species  which  are  significant  in  the  mass  and  pressure 
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balances  (46)  and  (47)  in  terms  of  and  the  less  significant 
species  in  terms  of  their  Jin  P^. 

The  Newton- Raphson  procedure,  as  applied  by  the  EST  pro¬ 
gram,  can  be  summarized  by  considering  a  set  of  equations  of  the 
general  form 


f  j  , ^2  9  •  •  •  ,  •  .  • )  —  0 

At  any  point  in  the  solution  procedure  there  exists  a  set  of  esti¬ 
mates,  x|,  for  all  the  variables  which  will  in  general  not  satisfy 
all  of  the  relations  and  will  lead  to  a  non-zero  value  of  the  f y 
namely,  e...  The  Newton-Raphson  method  proceeds  to  "drive"  these 
errors  toward  zero  by  evaluating  the  change  in  each  unknown  variable, 
Ax^,  which  would  reduce  all  the  errors  to  zero  if  the  functions,  f ^ , 
were  linear.  The  linear  approximation  is  based  on  the  current  values 
of  the  unknown  variables  and  the  corresponding  array  of  values  of 
the  partial  derivatives  Sfj/Bx^  Thus 


df 


j 


(50) 


which  is  locally  correct  and  is  integrated  to 


( Af  j)  * 


(51) 


in  the  linear  approximation.  The  solution  of  (50)  is 


dx. 

l 


(52) 


where  the  array  of  partial  derivatives  appearing  in  (52)  is  simply 
the  matrix  inverse  of  the  array  in  (50)  .  In  the  EST  program  the 
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formulation  of  the  partial  derivatives  uses  the  variables,  in  P^,* 
in  T  and  in#!,  and  (52)  yields,  for  example, 


d( in  Pi) 


3 

37 


in  P.. 

j 


(53) 


which  if  taken  as  linear  all  the  way  to  solution  yields 
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since  the  desired  change  in  the  functions  is  simply  the  negative  of 
the  error.  An  equally  exact  relation  obtained  from  (53)  is 
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which  if  taken  as  linear  all  the  way  to  solution  yields 


p  *  *  —  p  *  =  ^p  * 
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The  EST  program  uses  (56)  for  all  species  significant  in  mass  bal¬ 
ances  and  (54)  for  the  others. 


3.1.2  Restriction  on  Corrections 


The  set  of  correction  (Ax^)*  can  be  thought  of  as  a 
vector  in  the  space  of  the  independent  variables  which  is  added  to 
the  current  vector  approximation  xt  to  yield  a  new  estimate  x£*. 


* 

The  choice  of  in  P .  permits  a  matrix  reduction  by  the  use  of  the 
simple  algebraic  substitution,  previously  mentioned  after  (49), 
prior  to  matrix  inversion. 
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Experience  has  shown  that  it  is  frequently  unwise  to  proceed  along 
this  correction  vector  the  full  amount  indicated  by  (54)  or  (56) . 
Rather,  it  is  better  to  proceed  a  limited  way,  although  preserving 
the  same  direction.  At  other  times,  it  is  expedient  to  depart  from 
this  vector,  and  seek  another  based  on  freezing  the  value  of  some 
variable  and  eliminating  a  corresponding  equation. 

The  scaling  of  the  correction  vector  is  such  as  to  limit 
changes  in  the  partial  pressures  of  major  species  to  increases  of  one 
order  of  magnitude  and  decreases  of  three  orders  of  magnitude,  and 
changes  of  temperature  to  approximately  20  percent. 

Molecular  weight,  temperature  and  condensed  phase  concen¬ 
tration  corrections  are  frozen  and  a  new  correction  vector  generated 
if  the  initial  sat  of  corrections  indicate  excessive  temperature  or 
molecular  weight  excursions,  a  contradictory  temperature  change,  or 
negative  corrections  on  newly  introduced  condensed  species . 

The  formulation  of  these  and  other  scaling  and  freezing 
criteria  is  an  essential  feature  of  the  EST  program.  Because  of  these 
features,  convergence  is  virtually  assured  for  well  formulated,  physi¬ 
cally  unique  problems. 

3.1.3  Base  Species 

The  discussion  of  Section  2  described  the  equilibrium 
reaction  equations  as  equations  giving  the  formation  of  a  species 
from  atomic  elements.  Thus  the  reactants  are  elements  and  the  pro¬ 
ducts  are  usually  molecules.  This  scheme  has  the  advantage  of  for¬ 
mal  simplicity,  since  the  stoichiometric  coefficient  needed  in  the 
equilibrium  equations  are  given  directly  by  the  product  species  chem¬ 
ical  formula.  This  scheme  can  have  computational  disadvantages, 
however,  since  the  atomic  elements  are  frequently  not  present  to  any 
great  extent  in  the  equilibrium  system.  This  in  itself  results  in 
no  disadvantage.  If,  however,  a  molecule  (e.g.,  CO)  dominates  more 
than  one  mass  balance  (e.g.,  C  and  0)  loss  of  significant  figures  can 
slow  or  defeat  convergence. 
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It  is  more  desirable  to  write  the  equilibrium  reactions 
(1)  and  (2)  as  well  as  the  mass  balances  in  terms  of  reactant  species 
which  are  actually  present  in  appreciable  amounts.  These  species  are 
termed  "base  species"  (from  Ref.  20)  since  all  other  species  are  taken 
to  be  formed  from  them. 

The  EST  program  selects  the  base  species  from  the  candi¬ 
date  species  thermochemical  data  which  are  input  as  specified  in  Sec¬ 
tion  4.  The  program  automatically  selects  as  base  species  the  first 
set  of  species  satisfying  the  requirement  that  (1)  all  other  species 
may  be  formed  from  this  base  species  set  and  (2)  that  no  balanced 
reaction  can  be  written  involving  only  base  species.  One  base  spe¬ 
cies  may  be  considered  to  represent  each  element.  Thus,  the  base 
species  are  established  by  the  order  in  which  the  user  inputs  the 
candidate  species  thermochemicai  data.  The  calculation  of  the  stoi¬ 
chiometric  coefficients  and  equilibrium  constants  appropriate  to  any 
set  of  base  species  is  handled  automatically  by  the  program. 

While  the  selection  of  base  species  is  not  usually  cri¬ 
tical,  certain  computational  economies  may  be  realized  if  the  base 
species  are  those  species  which  occur  in  relatively  high  concentra¬ 
tions  in  the  system.  Thus,  the  user  is  advised  to  place  at  the  head 
of  the  species  thermochemistry  data  card  set  (Section  4)  a  group  of 
species  which  are  expected  to  occur  in  high  concentrations. 
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3.2  FORTRAN  ROUTINES 

In  order  to  permit  an  understanding  of  the  actual  solu¬ 
tion  mechanics  and  their  re] ation  to  the  FORTRAN  source  program,  a 
brief  description  of  each  of  the  subroutines  are  provided  in 
this  section.  For  convenience  the  routines  have  been  grouped, 
where  possible,  under  more  general  headings.  The  calling  of  these 
routines  and  the  principal  iterative  Loop  is  controlled  by  the  main¬ 
line  program.  All  routines  except  SQUEE  and  ETIMEF  are  called  by 
the  mainline  program. 

3.2.1  Thermodynamic  Data  Input  Routines 

IMELM:  Reads  elemental  data  information  and  normalizes 

elemental  composition  information. 

INPUT;  Reads  species  thermochemical  data,  selects  base 
species,  determines  stoichiometric  coefficients 
of  formation  reactions  of  all  other  species  from 
the  base  species ,  flags  condensed  species ,  sets 
up  first  guesses,  reads  diffusion  factor  and 
fail  temperature  data  (KR(3)^9)*.  and  calculates 
species  molecular  weights. 

BELCH:  Establishes  a  base  species-element  correspon¬ 

dence  table  which  is  required  if  an  input 
element  is  temporarily  omitted  from  the  chem- 
icaJ  system  being  analyzed. 

3.2.2  Problem  Setup  and  Initiation 

ZIPIN:  Reads  control  card,  relative  mass  or  ablation 

rate  data,  and  (if  requested)  title  information; 
controls  logic  for  temperature  fill-in  opera¬ 
tions  (KR (10) ) ;  and  initializes  and  sets  various 
control  flags. 

ft 

KR  (I )  refers  to  elements  of  an  input  array  which  controls  most  of 
the  program  options  -  see  Section  4. 
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ALPST :  Reads  diffusion  factor  values  if  KR(3)  =  7, 

8  or  9;  sets  diffusion  factors  based  on  expo¬ 
nential  approximation  or  input  values  using 
information  from  INPUT  and  ALPST;  determines 
which  elements  are  absent  from  the  system  and 
flags  corresponding  molecular  species  and,  in 
particular,  corresponding  base  species  (see 
BELCH) ;  reinitializes  species  omitted  from  the 
prior  solution  and  zeros  those  to  be  omitted 
from  current  solution;  evaluates  conserved 
quantity  parameters  entering  mass  balances . 

Calculation  of  Errors  and  Error  Derivatives  Used  Within 

Newton  Raphson  Iteration 

THERM:  On  first  iteration  and  after  every  change  of 

system  temperature  this  routine  is  called  to 
evaluate  molecular  thermodynamic  properties 
(e.g.,  enthalpy,  entropy,  free  energy,  specific 
heat)  and  the  equilibrium  constants  appropriate 
to  each  formation  reaction.  On  the  first  iter¬ 
ation  certain  reinitializations  are  performed 
including  assigning  temperature  to  3000 °K  if  the 
nonconvergent  flag  has  been  set  and  certain  key 
summations  are  evaluated.  Subsequently  these 
summations  are  determined  in  CRECT.  (The  rea¬ 
son  for  these  functions  inclusion  in  THERM 
relates  to  the  overlaying  of  subroutines  and 
molecular  species  data  in  an  IBM  1130  version 
of  this  program.) 

MATl:  Initializes  mass  balance  error  equations  and 

determines  contribution  of  base  species  to 
these  errors.  Initializes  matrix  of  partial 
derivatives  of  mass  balance  errors  with  respect 
to  log  of  base  species  partial  pressures  and 


o 
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introduces  base  species  contributions  to  these 
coefficients .  Commences  search  for  surface 
species .  On  first  iteration  normalizes  gas 
phase  partial  pressures. 

MAT2:  Evaluates  errors  in  formation  reaction  equilib¬ 

rium  relation  for  non-base  species.  Introduces 
contribution  of  non-base  species  to  mass  balance 
equations .  For  gases  the  contribution  of  the 
partial  derivatives  of  the  mass  balance  errors 
with  respect  to  the  non-base  species  log  par¬ 
tial  pressure  is  related  to  base  species  log 
partial  pressure  via  the  equilibrium  equation 
and  introduced  into  the  array  initialized  in 
MATl.  For  each  condensed  species  which  has 
been  previously  introduced  into  the  system  an 
additional  equilibrium  relation  is  added  to  the 
mass  balance  equations .  One  additional  con¬ 
densed  species  relation  is  accepted  if  equilib¬ 
rium  indicates  a  positive  formation  potential 
(see  AFMAT  relative  to  this  species)  .  Completes 
search  for  surface  species. 

MAT3:  Completes  formulation  of  mass  balance  errors  and 

their  derivatives .  Checks  for  convergence  of 
mass  and  equilibrium  equations. 

3.2.4  Calculation  of  Variable  Corrections 

Upon  inversion  of  the  matrix  of  error  derivatives  and 
its  multiplication  of  the  errors  (see  RERAY)  a  set  of  corrections 
are  obtained  for  temperature,  pressure-molecular  weight  product, 
log  partial  pressure  of  the  base  species  and  the  relative  moles  of 
condensed  species.  Before  these  constraints  can  be  applied,  how¬ 
ever,  a  significant  amount  of  messaging  is  required. 
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AFMAT:  Immediately  after  the  inversion,  performs  a 

series  of  tests  on  the  corrections  and  if 
necessary  imposes  constraints  on  equations  and 
repeats  inversion.  Examples  includes  the  newly 
introduced  condensed  species  with  zero  concentra¬ 
tion  has  a  calculated  negative  correction  -  the 
corresponding  equilibrium  equation  is  removed 
and  the  correction  set  to  zero;  if  the  tempera¬ 
ture  is  currently  equal  to  a  minimum  or  maximum 
value  dictated  by  phase  change  or  other  discon¬ 
tinuous  phenomena  and  if  a  contrary  temperature 
correction  is  predicted  -  the  temperature  con¬ 
straint  is  removed  and  the  temperature  frozen  - 
if  an  excessive  negative  change  in  the  pressure- 
molecular  weight  product  is  predicted  -  the 
pressure  constraint  is  removed  and  pressure 
frozen.  Certain  temperature  minima  and  maxima 
are  also  set  to  control  the  convergence  on  sur¬ 
face  equilibriu.ii  calculations. 

SCALE:  From  basic  corrections,  corrections  to  log  partial 

pressures  of  non-base  species  are  calculated. 

Scans  all  corrections  and  determines  maximum 
damping  factor  which  will  permit  all  corrections 
to  satisfy  certain  constraints.  These  include 
one  order  of  magnitude  increase  and  three  orders 
of  magnitude  decrease  in  partial  pressures  of 
significant  species,  and  less  restrictive  con¬ 
straints  on  less  important  species.  The  scale 
factor  generated  by  this  routine  will  subsequently 
be  applied  to  all  corrections.  Performs  correc¬ 
tion  of  temperature  and  pressure-molecular  weight 
product. 

CRSCT:  Performs  all  corrections  according  to  scale  fac¬ 

tors  calculated  in  SCALE.  Makes  corrections 
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based  either  on  linearization  of  mass  balances 
(delta  of  mole  fractions)  or  of  equilibrium  equa 
tions  (delta  of  log  of  mole  fractions)  depending 
on  relative  importance  in  mass  balances.  Evalu¬ 
ates  set  of  summations  previously  discussed  rela 
tive  to  THERM  in  Section  3.2.3. 

3.2.5  Calculation  of  Output  Quantities 

PROPS :  Determines  certain  derivative  properties  from 

the  inverse  matrix  of  error  derivatives  includ¬ 
ing  the  equilibrium  specific  heat  and  the  isen- 
tropic  exponent;  calculates  quantities  using 
diffusion  factor  model  and  outputs  calculated 
quantities . 

OUTPUT:  Generates  and  outputs  all  additional  terms  re¬ 
quired  for  output  as  displayed  in  Section  5, 

3.2.6  Utility  Routines 

RERAY:  Linear  equation  solver  and  inversion  routine. 

Call  list  is 

N  =  number  of  rows  in  C 
C  =  coefficient  array 
NN  =  number  of  columns  in  C  minus  N 
D  =  set  of  column  vectors  of  length  U 
NNN  =  number  of  column  vectors 
LS  =  column  rearrangement  flag  (zero  implies 
no  rearrangement) 

IS  =  Flag,  -2  at  call  results  in  before  and 
after  display  of  arrays,  <0  at  return 
implies  singularity  encountered  at  3-ow-lS. 
Returns  with  left  N  x  N  square  of  C  inverted 
and  thi3  inverse  multiplied  by  remaining 
columns  of  C  (if  any)  and  D. 
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SWAP; 


SQUEE: 


ETIMEF 

and 

ETIME: 


Replaces  one  array  and  vector  with  another  set, 
call  list  is 

A  array  1  (16  x  16  in  current  listing) 

B  vector  1  (16  in  length  currently) 

SLA  array  2 
SLB  vector  2 

Returns  with  A  and  B  replaced  by  SLA  and  SLB, 
respectively.  SLA  and  SLB  are  unchanged. 

Prepares  integer  and  alphanumeric  information 
from  real  variables  for  compressed  card  output. 
Call  list  is 

A  real  variable  to  be  converted 
I  first  of  three  storage  locations  for 
results 

IDEE  presumed  decimal  location  in  result 

Returns  with  1(1)  an  integer,  1(2)  the  sign  of 
the  exponent  (alphanumeric),  1(3)  the  exponent, 
such  that 


A 


Id) 


10 


IDEE 


10 


±1(3) 


System  time  routines  available  on  Univac  1108 
systems  tape.  For  other  systems  a  conversion 
or  dummy  routine  is  required.  Call  of  ETIME 
sets  zero  time  in  clock.  Subsequent  call  of 
ETIMEF (X)  retunrs  time  in  seconds  from  last 
call  of  ETIME  as  real  variable  X. 


3.3  DEFINITIONS  OF  SEVERAL  OUTPUT  QUANTITIES 

Several  quantities  appear  on  the  standard  program 
output  under  abbreviated  names.  In  this  section,  these  variables 
are  more  accurately  defined.  The  derivative  property  output  gives 
the  following: 

I  3h . 

CP-FROZEN  =  >  the  frozen  constant 

i=l  1  3T 

pressure  specific  heat  capacity 
(Btu/lbm°R  or  Cal/gm°K) 


CP-EQUIL  = 


ah\  I  3k.  v  I  3h.  . 

i)  -  LA-w  ■ 

p  1=1  p  1=1  'p 

the  equilibrium  constant  pressure 
specific  heat  capacity 
(Btu/lbm°R  or  Cal/gm°K) 


DLNM/DLNT  =  fgf) 

'  X ' 


(unitless) 


DLNM/DLNP  = 


(unitless) 


GAMMA  =  ^  -  Cp/Cv  i  t°r  an  ideal  gas 


(unitless) 


Other  property  output  includes : 


u1  =  £  x-iF-i '  a  9as  mixture 

j  3  J 

property  which  reduces  to  unity 
for  assumed  equal  diffusion  coef¬ 
ficients  Fj  =  F  =  1.0  (unitless) 
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MU 2  *»  y2  =^^!jXj/F j  ,  a  gas  mixture 
property  which  reduces  to 
the  mixture  molecular 
weight, ^  ,  for  assumed  equal 
diffusion  coefficient 
F j  =  F  =  1.0  (grams  of  sys¬ 
tem/mole  of  gas) 

M0L.WT  =  W!  ,  the  composite  mixture 

molecular  weight  (see  Equation 
6,  page  2-4)  ( grams  system/moles 
of  gas) 

I 

HT1L  =  fi  =  £  Z . h . ,  a  property  of  the 

j-1 

gas  mixture  which  reduces  to  the 
static  enthalpy,  h,  for  assumed 
equal  diffusion  coefficients. 

Z^  is  defined  by  Equation  36, 
page  2-21.  (Btu/lbm) 

I  3h. 

CPTIL  =  C  =  £  Z,-5m  ,  a  property  of 

p  i=l  p 

the  gas  mixture  which  reduces  to 
the  frozen  specific  heat  capacity 
for  assumed  equal  diffusion 
coefficients.  is  defined  by 
Equation  36,  page  2-21 
(Btu/lbm°R  or  Cal/gm°K) 

I 

HTIL*  =  h*  =  X]  Z?h.,  a  property  of  the 
i=l  1  1 

gas  mixture  which  reduces  to 
HTIL  when  GAMEX  (the  unequal 
diffusion  exponent)  -  1.0,  and 
reduces  to  static  enthalpy,  h, 
for  assumed  equal  diffusion  co¬ 
efficients.  Z*  is  defined  by 
Equation  35,  page  2-21.  (Btu/lbm) 


3-15 


GAMEX  =  y,  the  weight  factor  exponent 
between  K  and  Z  mass  fractions 
used  to  obtain  the  Z|  values 
in  Equation  35,  page  2-21. 
(unitless) 

ELEMENTAL  K  MASS  FRACTIONS  =  Kk,  usual  definition  of  mass 

fraction  for  element  K  ir»  the 
system,  see  Equation  10,  page 
2-10 .  (gm  of  k/gm  system) . 

ELEMENTAL  Z  MASS  FRACTIONS  =  Z*  ,  elemental  fT  weighted  mass 

fraction  of  element  K  in  the 
system,  defined  by  Equation  34, 
page  2-21.  (gm  of  K/gm  system) . 

All  of  the  above  properties  are  calculated  and  output 
for  each  problem  solution  using  the  current  values  of  I\  and  GAMEX, 
even  though  the  particular  option  being  computed  does  not  require 
these  values .  An  important  example  is  the  equal  diffusion  surface 
equilibrium  option.  Even  though  the  HTIL*  value  is  computed  and 
printed  for  a  GAMEX  value  of  0.667,  the  value  of  HTIL*  punched  on 
output  cards  (for  input  to  a  heat  conduction  program)  will  be 
static  enthalpy,  h,  as  it  should  be. 

Additional  descriptions  of  the  variables  computed  and 
output  by  the  EST  code  are  given  in  References  2,  3,  and  4. 

3.4  PROGRAM  OPERATION  DETAILS 

3.4.1  Program  Storage  Requirements 

The  current  version  of  the  EST  program  allows  149  chem¬ 
ical  species,  10  elements,  a  total  of  up  to  16  simultaneous  equations 
(1  pressure  balance,  10  mass  balances,  4  present  condensed  species, 

and  1  enthalpy,  temperature  or  surface  equilibrium  condition),  50  B' 

y 

and  B'  entries,  and  50  entries  in  the  diffusion  factor  tables.  With 
c 

these  limits  approximately  25,000  (decimal)  words  of  storage  are 
required. 


3-16 


3.4.2  Tape  Requirements 

All  input  and  output  tape  designations  are  set  at  the 
beginning  of  the  main  program  with  one  exception.  The  designations 
are 


KIN 

System  input 

(=5) 

KOUT 

System  output 

(=6) 

JAN 

Scratch  (or  save) 

(=18) 

tape,  see  KR(3) 

In  RERAY,  KOUT  is  also  set  to  the  system  output  tape 
(=6  currently) . 

3.4.3  Operator  Controls 

No  operator  intervention  control  switches  exist.  All 
controls  are  set  at  the  initia ;ion  of  execution  of  each  problem  with 
the  corresponding  control  card  (see  Section  4) . 

3.4.4  Dumps  and  Diagnostics 

Control  column  7  (KR{7))  controls  the  output  of  re¬ 
quested  diagnostic  data.  The  interpretation  of  this  data  is  not 
always  obvious  and  users  should  consult  with  the  program  authors 
if  interpretation  is  mandatory.  This  same  data  will  automatically 
be  obtained  if  more  than  67  iterations  are  required  for  conver¬ 
gence  and  will  continue  through  the  70th  and  last  iteration.  Other 
diagnostic  complaints  may  be  self-explanatory.  The  failure  to 
obtain  convergence  is  usually  related  to  physical  situations  for 
which  no  unique  solution  exists.  For  example  specification  of  a 
graphite  B*  less  than  about  0.09  in  equilibrium  with  air  results 
in  the  program  flailing  about  the  program  minimum  temperature  (300°K) 
since  the  C02  plateau  represents  the  minimum  equilibrium  B'. 

3.4.5  Program  Running  Time 

Program  run  time  is  most  conveniently  measured  in  terms 
of  the  number  of  solutions,  and  time  per  solution,  required  for  a 


a 
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given  problem.  Time  per  solution  is  conveniently  measured  in  terms 
of  the  number  of  iterations,  and  time  per  iteration,  required  for  a 
converged  solution.  Thus,  the  run  time  for  a  given  problem  may  be 
estimated  as : 


Run  Time  =  No. 


of  Solutions  x 


Iterations  Time 
Solutions  x  Iteration 


The  number  of  solutions  required  for  a  given  problem  is 
determined  by  the  user  (except  when  one  of  the  fill-in  options  are 
used,  in  which  case  the  exact  number  of  inserted  or  appended  solu¬ 
tions  is  not  known  a  priori) .  For  example,  for  simple  closed  sys¬ 
tem  problems,  the  number  of  solutions  is  equal  to  the  number  of 
times  card  set  1  (Section  4)  is  repeated,  and  for  surface  equilib¬ 
rium  open  system  problems,  the  number  of  solutions  is  equal  to  the 

number  of  elements  in  the  input  B*  x  B'  matrix. 

r  g  c 

The  number  of  iterations  required  for  a  converged  solu¬ 
tion  is,  of  course,  a  function  of  the  type  of  problem  being  solved. 
The  sample  problems ,  Section  5 ,  should  provide  some  guidance  in 
estimating  this  quantity.  Roughly,  most  solutions,  other  than  the 
initial  solution,  converge  in  about  5  to  15  iterations. 

The  time  required  per  iteration  is  a  function  of  the 
computing  machine  as  well  as  the  problem  type.  Machines  in  the 
speed  range  of  the  Uni vac  1108  or  CDC  6600  require  on  the  order  of 
10  to  60  milliseconds  per  iteration.  Again,  the  sample  problems 
in  Section  4  (which  were  run  on  a  Univac  1108)  should  provide  some 
guidance  here. 
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SECTION  4 
INPUT  PREPARATION 

This  section  defines  the  format  and  significance  of  each 
field  of  the  input  data  card  deck  for  Version  3  of  the  Aerotherm 
Equilibrium  Surface  Thermochemistry  (EST)  program.  The  input  con¬ 
sists  of  six  card  sets.  These  are: 

1 .  CONTROL  CARD 

2 .  TITLE  INFORMATION 

3 .  SURFACE  FLUXES 

4.  ELEMENTAL  COMPOSITION  DATA 

5.  DIFFUSION  FACTOR  DATA 

6.  SPECIES  THERMOCHEMICAL  EQUILIBRIUM  DATA 

Card  Set  1  is  a  control  card  which  contains  a  iO-element 
array,  KR(I) .  This  array  controls  most  of  the  program  options  and 
tells  the  program  what  to  expect  from  the  remaining  card  sets .  Few 
problems  require  all  card  sets.  The  card  sets  are  described  in  de¬ 
tail  in  this  section. 
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CARD  SET  1  -  CONTROL  CARD  (1  Card) 

FIELD  1  (Columns  1-10,  Format  1011)  this  is  the  variable 
array  KR(I)  which  is  used  to  control  the  various 
program  options . 

COLUMN  1  —  STATE  OPTION 

0  Assigned  temperature 

1  Surface  equilibrium 

2  Assigned  enthalpy 

3  Assigned  temperature  supressing 

surface  equilibrium 
8  Stop  Program  Execution 

COLUMN  2  —  ELEMENTAL  COMPOSITION 


0  Accept  the  resident  values 

(i.e.,  those  used  in  the 
previous  solution) 

1  Input  new  elemental  composi¬ 

tions  for  the  system  compo¬ 
nents  (card  set  4) 

COLUMN  3  —  SPECIES  THERMOCHEMICAL  DATA 


0  9 


1234 


34 


5678 

789 
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Accept  the  resident  values  for 
species  thermochemistry  data 
(i.e.,  those  used  in  the  pre¬ 
vious  solution) 

Read  from  cards  new  species 
thermochemistry  data  (card 
set  6)  and,  if  provided,  data 
on  diffusion  factors  (card 
set  5) 

Store  all  species  data  including 
diffusion  factors  read  above 
from  card  sets  5  and  6  on  tape 
Read  all  species  data  from  tape 
as  last  stored 

Update  values  of  diffusion  fac¬ 
tors  (card  set  5) 

Print  out  the  thermochemistry 
data  for  all  species  read  or 
stored  which  may  be  formed  from 
the  elements  input  (card  set  4) 


COLUMN  4  --  MASS  BALANCE  OPTIONS 


0  Closed  system  mass  balance 

1234  Open  system  mass  balances 

1  Equal  diffusion  coefficients.  This 

flag  overrides  the  use  of  all  un¬ 
equal  diffusion  coefficient  equations. 
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COLUMN  4 

234 

2 

3 

4 

COLUMN  5 
0 

1 

COLUMN  6 
COLUMN  7 
0 

1 

j 


COLUMN  8 

0 

1 

2 


—  MASS  BALANCE  OPTIONS  (Continued) 

Hence,  no  corrections  to  F. 

l 

vaues  or  GAMEX  are  required. 
Unequal  diffusion  coefficients 
(requires  boundary  layer  edge 
composition) 

No  change  from  previously  saved 
boundary  layer  edge  composition 
Use  the  prior  solution  as  the 
boundary  layer  edge  composition 
Determine  the  boundary  layer  edge 
composition  ,  in  terms  of  Z  mass 
fractions ,  from  the  elemental 
data  cards  (card  set  4)  as  last 
read. 

—  TITLE  INFORMATION  OPTION 

No  effect 

Read  three  cards  and  print  contents 
on  the  top  of  problem  output. 

—  Not  used 

—  DIAGNOSTIC  OUTPUT  CONTROL 

No  diagnostic  output 
Output  a  single  line  of  diagnostic 
information  per  iteration 
(where  j  is  a  digit  greater  than 
one)  Output  full  diagnostic  in¬ 
formation  for  5 ( j— 1)  iterations 
and  output  a  single  line  of  diag¬ 
nostic  information  for  subsequent 
iterations . 

—  PUNCHED  CARD  OUTPUT  FORMAT 
(implemented  only  if  KR(10)  ^  0) 

Format  appropriate  to  input  the  CMA 
program 

Format  appropriate  to  input  the 
ASTHMA  program 

Output  cards  in  both  of  the  above 
formats 


COLUMN  9  —  Not  used 
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COLUMN  10  —  OPEN  SYSTEM  SOLUTION  CONTROL  AND  PUNCHED 
CARD  OUTPUT  OPTIONS  (used  to  provide 
punched  card  output  for  input  to  energy 
balance  programs) 


0 

123456789 


2345 


6789 
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45  89 


No  effect 

For  each  convergent  solution,  out¬ 
put  a  card  or  cards  of  data  appro¬ 
priate  to  energy  balance  calcula¬ 
tions  (see  KR ( 8) ) 

Input  an  array  of  pyrolysis  gas 
(component  2)  and  char  (component 
3)  fluxes  (card  set  2)  .  For  each 
pyrolysis  gas  flux,  obtain  a  sol¬ 
ution  for  all  char  fluxes  (valid 
only  for  KR(1)  =  1) 

Using  last  input  arrays,  obtain  sol¬ 
utions  as  above  (valid  only  for 
KR(1)  =  1) 

For  each  pyrolysis  gas  rate,  append 
a  sequence  of  assigned  temperature 
solutions  (suppressing  surface 
equilibrium)  at  the  last  char  flux 
commencing  at  the  next  500 multi¬ 
ple  below  the  last  surface  equilib¬ 
rium  temperature  and,  after  a  suc¬ 
cession  of  500 °K  decrements,  termin¬ 
ate  at  500 °K 

For  each  pyrolysis  gas  rate,  insert 
up  to  ten  additional  char  fluxes 
between  each  pair  of  tabular  char 
flux  entries  as  required  to  limit 
temperature  gaps  between  solutions 
to  less  than  500 °K 


FIELD  2  (Columns  11-20,  Format  F10.5)  The  thermodynamic  state 
variable,  Z 


If  KR(1)  =  0,  1,  3,  and  Z  ^  0:  Z  is  assigned  to  temperature 

(°K) 

If  KR(1)  =  2  and  Z  /  0:  Z  is  assigned  to  enthalpy  (cal/gm) 

If  Z  =  0 :  resident  values  of  temperature  or  enthalpy  (which¬ 
ever  is  consistent  with  the  KR(1)  specification)  will  be  used 

A  negative  temperature  entry  results  in  a  frozen  composition 
of  the  system  (as  previously  calculated)  and  the  calculation 
of  properties  of  this  system  at  the  absolute  value  of  the 
assigned  temperature 

If  KR(1)  =*  1,  the  temperature  assignment  serves  only  as  a 
first  guess  supplanting  the  resident  temperature  (on  the 
first  solution  or  after  completing  solutions  for  an  array 
of  pyrolysis  gas  and  char  rates,  the  first  guess  tempera¬ 
ture  is  taken  as  3000 °K) 


o 
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FIELD  3  (Columns  21-30,  Format  F10 . 5)  Pressure  variable,  PR 

If  PR  ^  0,  PR  is  assigned  to  the  system  pressure  (atmospheres) 
If  PR  =  0,  the  resident  value  of  pressure  will  be  used 

Field  4  (columns  31-40,  Format  F10.5)  the  Z  mass  fraction  weight¬ 
ing  factor  exponent,  GAMER 

7  f  GAMER  /  0,  GAMER  is  assigned  to  GAMEX  which  is  the  weight¬ 
ing  factor  exponent  for  the  mass  transfer  driving  potential 
fractions  (Z*)  in  terms  of  the  mass  and  Z  fractions  where 

Z*  a  k1-GAMEX 

If  GAMER  =  0 ,  no  change  will  be  made  in  the  weighting  factor 
exponent  GAMEX.  GAMEX  is  set  equal  to  2/3  prior  to  the  first 
solution,  and  this  value  is  appropriate  for  most  calculations 

Field  5  (Columns  41-70,  Format  3F10.5)  Relative  component  mass 
variables  WR(1) ,  WR(2) ,  WR(3) 

If  any  of  these  three  variables  is  non-zero,  these  variables 
are  assigned  to  W(l) ,  W(2) ,  W(3) ,  the  relative  masses  of 
components  1,  2,  and  3  respectively. 

For  surface  mass  balances ,  the  relative  masses  are  the 

relative  mass  fluxes  normalized  by  p  u  C„  and: 

e  e  M 

W(l)  =  boundary  layer  edge  gas 
W (2)  =  pyrolysis  gas 
W  ( 3)  =  char 

If  KR(10)  >  1,  this  field  will  be  ignored. 

If  all  WR  are  zero,  resident  values  of  W  are  used  (prior  to 
the  first  solution,  W(l)  =  1.0,  W(2)  =  0,  and  W(3)  =  0). 

Field  6  (Columns  71-80,  Format  2A4 ,  A2)  Input  a  title,  job  name, 

etc.,  which  will  appear  on  all  printed  and  punched  output. 


CARD  SET  2  -  TITLE  INFORMATION  (3  cards) 

IF  KR(5)  =0  skip  this  card  set 

Cards  1,  2,  and  3  (Columns  1-80,  Format  20A4) 

Enter  any  80  alphanumeric  characters  on  each  card. 
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-  SURFACE  FLUXES 


(number  of  cards  =  3  + 
B^'s  -1)  +  1/8  x  (no. 


1/8  x  (no,  of 
of  B”s  -1) 


If  KR(10)  /  2,3,4  or  5,  skip  this  card  set 

If  KR(10)  =  2,3,4,  or  5,  the  normalized  pyrolysis  gas  and  char 
fluxes  are  input  here.  These  fluxes  are  normalized  with  respect 
to  the  mass  transfer  coefficient  as  per: 


pyrolysis  gas  surface  mass  flux 


P  uTT. 
He  e  M 


B 


g 


char  surface  mass  flux 
peueCM 


CARD  1 

FIELD  1  (Columns  1-3,  Format  13)  Number  of  B’’s  to  be  read 
from  the  next  card(s)  g 

FIELD  2  (Columns  4-6,  Format  13)  Number  of  B*'s  to  be  read 
from  subsequent  card(s) 

CARD (S)  2  (Columns  1-80,  Format  8F10.)  Enter  values  of  B'*s,  eight 
to  a  card,  to  a  maximum  of  50  g 

For  convenience  of  subsequent  energy  balance  input 
(using  card  output  from  this  program) ,  the  B' 's  should 
be  ordered  in  either  ascending  or  descendinggsequence 

CARD ( S )  3  (Columns  1-80,  Format  8F10.)  Enter  values  of  B^'s,  eight 
to  a  card 

These  values  should  be  in  descending  sequence 


NOTE :  Care  should  be  taken  to  avoid  zero  entries  if 
by  so  doing  elements  are  not  available  from 
which  condensable  species  may  be  formed.  If 
edge  gas  is  air,  simultaneous  zero  assignments 
for  pyrolysis  gas  and  char  flux  reduces  the 
the  system  to  nitrogen  and  oxygen;  i.e.,  no 
condensables.  Note  that  from  inert  pyrolysis 
gases,  no  condensables  are  available. 
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CARD  SET  4  -  ELEMENTAL  COMPOSITION  DATA  (Number  of  cards  =  number 
.  ^  of  elements  +  1) 

If  KR(2)  =  0,  skip  this  card  set 


CARD  1  (Columns  1-3,  Format  13)  The  number  of  elements  in  the 
system  (in  ionized  systems,  the  electron  is  counted  as 
an  additional  element) 


CARD  2,3,4  .  .  .  .  (one  such  card  for  each  element) 

FIELD  1  (Columns  1-3,  Format  13)  the  atomic  number  iden¬ 
tifier  of  the  element  (the  electron  which  must 
be  included  for  ionized  closed  system  calcula¬ 
tions  must  have  identifier  99.  Any  non-zero 
identifier  may  be  used  for  other  elements  in 
this  list  but  thermochemical  species  formed 
from  this  element  must  have  an  identical  iden¬ 
tifier) 

FIELD  2  (Columns  4-15,  Format  3A4)  the  name  of  the 
element  (for  output  identification  only) 

FIELD  3  (Columns  16-25,  Format  F10.)  the  atomic 
weight  of  the  element 

FIELD  4  (Columns  26-55,  Format  3F10.)  the  relative 
amounts  of  the  element  in  components  1,  2, 
and  3,  respectively  (see  field  5  of  card  set  1) 

Positive  values  are  in  relative  gram-atomic- 
weights  (or  moles) 

Negative  values  are  in  relative  masses 

If  KR(4)  =  4  (column  4  of  card  set  1) ,  columns 
26-35  should  vontain  the  Z  mass  fraction  of  the 
element  in  che  boundary  layer  edge  gas .  They 
must  be  entered  with  a  minus.  These  values  are 
most  conveniently  obtained  from  a  prior  solu¬ 
tion  where  they  are  output  corresponding  to  the 
relevant  atomic  numbers.  Identical  diffusion 
factors  must  be  used  to  insure  meaningful  results. 

If  elements  are  added  to,  deleted  from,  or  rear¬ 
ranged  with  respect  to  a  prior  input,  it  is  nec¬ 
essary  to  re-read  card  set  6  either  directly  or 
from  the  tape. 
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FIELD  4  (Continued) 

A  typical  example  of  Field  4  input  is  as 
follows:  If  component  3  consists  of  50% 
mass  of  C  and  50%  mass  of  SiC^,  or  equiva¬ 
lently,  60.06  gm  of  C  for  each  60.06  gm  of 
Si02  (1  mole  of  SiOj) ,  or  finally  60.06  gm 

of  C  for  each  1.  gram-atomic-weight  of  Si 
and  2  gram-atomic-weight  of  0 ,  it  would  be 
most  convenient  to  input  -60.06,  2.0,  and 
1.0  into  columns  of  36-45  of  the  cards  for 
the  elements  C,  0,  and  Si. 
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CARD  SET  5  -  DIFFUSION  FACTOR  DATA  (number  of  cards  =  1/4  x  (7  +  no. 
- — — —  of  data  items,  viz.,  diffusion  factors,  reference  molec¬ 
ular  weight,  diffusion  exponent) 


If  KR(3)  =  5  or  6 ,  this  card  set  must  be  skipped 
If  KR( 3)  =  1,2,3  or  4,  this  card  set  may  be  skipped 
If  KR( 3)  =  7,8  or  9,  this  card  set  must  not  be  skipped 

CARD  1  (Columns  1-3,  Format  13)  the  total  number  of  data  items 
(viz.,  diffusion  factors,  reference  molecular  weight, 
diffusion  exponent)  to  be  entered 

Diffusion  factors  may  be  specified  for  any  or  all 
species  individually  and/or  diffusion  factors  may  be 
calculated  via  the 

.  /Molec.Wt.  \FFA 
Diffusion  Factor  =  - ) 

where  REFM  and  FFA  may  be  specified  here.  For  species 
for  which  diffusion  factors  are  not  specified  individ¬ 
ually,  diffusion  factors  will  be  calculated  via  the 
above  correlation,  and  if  REFM  and  FFA  are  not  speci¬ 
fied,  diffusion  factors  will  be  calculated  from  the  above 
with  REFM  =  23.4  and  FFA  =  0.431,  assigned  internally. 

CARD ( S )  2 


FIELDS  1, 3,5,7  (Columns  1-8,  21-28,  41-48,  61-68,  each 
format  2A4) 

The  chemical  symbol  of  the  species  for 
which  data  is  to  be  provided,  exactly  as 
it  appears  in  columns  73-80  of  the  first 
card  of  the  3  thermochemical  data  cards 
(card  set  6)  for  that  species.  To  spec¬ 
ify  values  of  REFM  or  FFA  in  the  next 
field,  enter  the  alpha  characters  "REFM" 
of  "FFA"  respectively  here. 

FIELDS  2, 4, 6, 8  (Columns  9-20,  29-40,  49-60,  69-80,  each 
format  E12 . 4) 

If  the  name  of  a  chemical  species  was 
entered  in  the  preceding  field:  the 
number  in  this  field  is  presumed  to 
be  a  diffusion  factor. 

If  the  name  REFM  or  FFA  was  entered  in 
the  preceding  field:  enter  the  desired 
values  of  REFM  or  FFA  respectively  in 
this  field. 
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The  diffusion  factors  introduced  into  the  solution 
in  this  fashion  supplant  the  values  set  automatically 
via  the  above  correlation,  with  REFM  =  23.4  and  FFA  = 
0.431,  respectively)  every  time  species  thermochemical 
data  are  read  from  cards  or  tape  (i.e.,  for  KR(3)  be¬ 
tween  i  and  8  inclusive) . 


This  card  set  data  is  not  saved  in  core.  However, 
this  card  set  will  be  stored  on  tape  together  with 
the  current  card  set  6  if  KR(3)  =  3  or  4,  and  thus  will 
be  introduced  automatically  if  the  tape  is  subsequently 
read  with  KR(3)  =  5  or  6.  If  KR(3)  =  7  or  8,  an  addi¬ 
tional  card  set  5  is  expected  and  will  be  temporarily 
appended  to  any  card  set  5  stored  on  tape.  If  this 
composite  set  contains  competing  data  for  the  same 
species,  the  latter  entry  will  be  chosen.  Latter 
entries  of  REFM  and  FFA  will  be  chosen.  If  KR(3)  = 

9,  only  those  species  indicated  on  the  card  set  will 
be  modified. 


o 
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CARD  SET  6  -  SPECIES  THERMOCHEMICAL  EQUILIBRIUM  DATA  (Number  of 
cards  =  1  +  3  x  number  of  species) 

If  KR(3)  =  0,5, 6, 7,8  or  9,  skip  this  card  set 

There  are  three  of  these  cards  for  each  molecular,  atomic,  condensed, 
ionic,  or  electron  species.  The  end  of  this  card  set  is  signaled  by 
a  blank  card  (see  reference  3  for  additional  discussion  relative  to 
this  card  set) 


CARDS  1,4,7  .  .  .  Describe  the  elemental  composition  of  the 

species  and  establish  its  name  designation 


FIELDS  1.3,5, ... ,13.  (one  for  each  element  in  the  species) 

(Columns  1-3 , 7-9 , 13-15 ,..., 37-39  , 
each  format  F3.0) 

The  number  of  atoms  (of  atomic  number 
given  in  the  following  field)  in  a 
molecule  of  this  species 

(If  field  1  is  zero,  this  card  is  pre¬ 
sumed  to  represent  the  end  of  card 
set  6) 


FIELDS  2,4,6, ... ,14.  (one  for  each  element  in  the  species) 

(Columns  4-6,10-12,16-18, . . . ,40-42, 
each  format  13) 

The  atomic  number  identifiers  (cor¬ 
responding  to  Card  Set  4)  of  the  ele¬ 
ments  in  the  molecule  (the  number  of 
atoms  of  which  was  given  in  the  pre¬ 
vious  field) 

(If  field  2  is  zero,  this  card  is  pre¬ 
sumed  to  be  the  first  card  of  card  set  5) 

Ionized  species  are  described  by  the  addition  or  sub¬ 
traction  of  an  electron  (atomic  number  99).  For  example, 

NO+  would  be  described  in  fields  1  through  6  as 
bb lbb 7bblbb8b-lb99  . 


FIELD  15  (Columns  43-72,  Format  7A4A2)  the  source  and  date 
of  the  thermochemical  data  for  this  specie.  Used 
for  output  only 


FIELD  16  (Columns  73-80,  Format  2A4)  the  name  designation  of 
this  species  (e.g.,  AL203) .  This  variable  is  used 
for  output  and  as  a  means  of  identifying  data  en¬ 
tered  in  card  set  5.  Typically  the  species  chemical 
symbol 
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CARDS  2 ,5, 8. . .Lower  temperature  range  thermochemical  data 


FIELDS  1,2, 3, 4 ,5  and  6  (Columns  1-54,  Format  6E9.6)  Input 

the  5  constants  (FI ,F2 ,F3 ,F4 ,F5 ,F6) 
appropriate  to  the  lower  temperature 
range  of  the  thermodynamic  data  for 
this  species.  These  constants  are 
defined  as  follows ,  where  T  is  in 
°K: 

FI  =  the  heat  of  formation  of  the 

species  at  298°K  from  the  JANAF 
base  state  (elements  in  most 
natural  form  at  298°K)  in  cal/ 
mole 

F2  =  the  enthalpy  change  of  the 

species  from  298°K  to  3000°K  in 
cal/mole 

F3,F4  and  F5  are  defined  by  a  curve 
fit  to  the  heat  capacity  at  constant 
pressure  of  the  form: 

Cp  =  F3  +  F4  T  +  p-  in  cal/mole  °K 

F6  =  the  entropy  of  the  species  at 
3000°K  in  cal/mole  °K 


FIELD  7  (Columns  61-66,  Format  F6.0)  the  upper  limit  of  the 
lower  temperature  range  in  °K 


FIELD  8  (Column  67,  Format  11)  the  phase  specification: 

1  signifies  gaseous  species 

2  signifies  solid  species 

3  signifies  liquid  species 


The  only  phase  combination  allowed  in  one  three  card  set 
is  solid-liquid  in  which  case  a  2  and  3  would  appear  in 
fields  8  of  cards  2  and  3,  respectively. 


o 
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CARDS  3,6 ,9 .Upper  temperature  range  thermochemical  data 

These  cards  are  the  same  as  cards  2,5,8...  except 
use  constants  for  the  upper  temperature  range  and 
field  7  is  ignored 


LAST  CARD  OF  CARD  SET  6  MUST  BE  BLANK 

The  end  of  the  species  thermochemical  data  is  signified  by  a 
blank  card.  Hence,  the  last  card  of  card  set  6  must  be  a 
blank  card. 


The  arrangement  of  these  three  card  sets  is,  mathematically, 
unimportant  for  equilibrium  solutions.  In  fact  if  no  condensed 
phase  species  exist  within  the  system  the  convergence  path  should 
be  identical.  Numerically  this  is  not  always  true.  A  single 
species  can  dominate  more  than  one  mass  equation,  which  in  turn 
can  lead  to  singular  arrays.  This  can  be  avoided  by  placing  these 
dominate  species  at  the  front  of  the  deck. 

The  program  automatically  selects  as  "base  species"  the  first 
set  of  species  satisfying  the  requirement  that  (1)  all  other 
species  may  be  formed  from  this  base  species  set  and  (2)  that  no 
balanced  reaction  can  be  written  involving  only  base  species . 

One  base  species  may  be  considered  to  represent  each  element.  In 
the  numerics  of  the  program,  the  base  species  actually  become 
the  "elements"  of  the  mass  balances.  The  means  of  establishing 
whether  one  species  may  be  formed  from  another  are  the  specifica¬ 
tions  given  in  fields  1  to  14  of  cards  1,4,7,....  Any  reaction 
which  can  occur  based  on  these  specifications  will  be  considered 
equilibriated  in  the  system  under  study.  To  suppress  equilibrium 
it  is  thus  required  to  change  names  (i.e.,  atomic  numbers)  of. 
certain  elements  and  to  introduce  new  "elements." 
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STOP  OR  CONTINUE 


The  program  normally  expects  problems  or  cases  to  be  stacked. 

An  8  punch  in  column  1  signifies  the  end  of  an  input  data  deck. 
Thus,  if  the  program  is  to  exit  in  a  graceful  fashion,  the  last 
card  in  the  input  data  deck  should  have  an  8  puncn  in  column  1. 
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SECTION  5 
SAMPLE  PROBLEMS 

Presented  in  this  section  are  five  sample  problems 
which  were  run  on  a  Univac  1108.  These  problems  are  relatively 
simple  but  do  exercise  most  of  the  EST  program  options.  The 
table  on  the  following  page  gives  a  brief  discussion  of  the 
nature  of  each  sample  problem.  In  addition,  for  each  sample 
problem,  the  following  is  presented: 

•  A  listing  of  the  input  data  deck 

•  A  few  typical  pages  of  the  program  output 
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DESCRIPTION  OF  EST 
SAMPLE  PROBLEMS 


Description 

No.  1  Closed  system  calculations  of  the 
thermochemical  states  of  two 
octane-air  ratios  ("Stoich? ometric" 
and  "llu2  theoretical  air) ,  at  one 
pressure  (10  atm)  and  four  temper¬ 
atures  (3000,  2500,  2000,  500°K) . 
119  different  molecular  species  are 
considered,  including  condensed 
carbon  (C*j . 

No.  2  Closed  system  calculation  of  the 
thermochemical  states  of  an  alu¬ 
minized  solid  rocket  propellant. 

The  propellant  consists  of  70.4 
percent  ammonium  perchlorate  (com¬ 
ponent  1),  9.6  percent  butadiene 
binder  (component  2)  and  20  percent 
aluminum  (component  3) .  Specified 
enthalpy  calculations  were  per¬ 
formed  at  two  conditions  (-777.4 
cal/gm,  68  atm  and  -1369.1  cal/gm, 

2 . 0  atm) . 


Sample 

Problem 


No.  3  Closed  system  calculation  of  the 
thermochemical  states  of  air  at 
one  pressure  (1.0  atm)  and  five 
temperatures  (2000,4000,6000, 

10 , 000°K) .  Elemental  composition 
includes  nitrogen,  oxygen,  argon, 
and  the  electron. 

No.  4  Open  system  calculation  of  the 

equilibrium  surface  thermochemi¬ 
cal  response  of  carbon  exposed 
to  air.  Equal  diffusion  coeffi¬ 
cients  are  employed  and  an  array 
of  B^,har  values  (2.0,  1.0,  0.5, 

0.2,  0.18,  and  0.17)  is  input. 

The  same  B^ar  values  are  utilized 
for  three  pressures  (1.0,  10.0  and 
100.0  atm)  and  the  option  to  in¬ 
sert  up  to  10  B^ar  values  is 
employed. 

No.  5  Closed  and  open  system  calculations 
of  the  surface  response  of  carbon 
(or  graphite)  phenolic  to  air.  An 
assigned  temperature  closed  system 
edge  calculation  (5530°K,  0.1  atm) 
defines  the  edge  gas  composition, 
a  frozen  edge  gas  table  is  gener¬ 
ated  (7  temperatures) ,  and  two  sets 
of  matrices  are  input. 


This  problem  demonstrates  the  use 
of  two  component  input  for  a  closed 
system  calculation.  Air  is  input 
as  component  1  (mass  fractions  of 
oxygen  and  nitrogen) ,  and  octane 
(c8h18)  input  as  component  2. 
Relative  mass  of  each  component  is 
input  on  card  set  1.  Note  that 
blanks  imply  reuse  of  previous  data. 

Three  components  are  used  to  define 
the  elemental  quantities  in  this 
problem.  Condensed  phase  alumina 
(AI2O3*)  is  found  to  be  present  in 
the  system  and  the  resultant  output 
format  is  shown.  For  the  -1369.1 
cal/gm,  2.0  atm  calculation,  the 
program  determined  that  25.452  per¬ 
cent  of  the  condensed  AI2O3*  was 
liquid.  The  option  to  read  and 
print  out  three  )ines  of  title  in¬ 
formation  is  also  demonstrated  in 
this  problem. 

This  problem  demonstrates  the  com¬ 
putation  of  the  thermochemical 
states  of  disassociated  and  ionized 
species.  The  output  of  species 
thermochemical  data  is  also  shown  in 
this  sample  problem. 


This  problem  demonstrates  the  calcu¬ 
lation  of  an  equal  diffusion,  equil- 
brium  surface  thermochemistry  table. 
Several  B^.har  values  were  inserted 
between  the  specified  values  to  limit 
the  maximum  temperature  gap  in  the 
table  tr  less  than  500°K.  The  typi¬ 
cal  output  is  from  the  10  atm  table. 
Also  shown  is  the  CMA  format  punched 
card  output  (also  written  on  printed 
output) . 


This  problem  shows  the  typical  tech¬ 
niques  used  to  generate  a  matrix  of 
Bgas  ani3  Bchar  values  suitable  for 
input  to  the  CMA  code.  The  effects 
of  unequal  diffusion  coefficients 
(defined  internally)  are  accounted 
for  and  nonablating  table  entries  are 
appended  to  the  end  of  the  final  set 
of  B£har  values •  Note  that  this  does 
not  represent  a  complete  surface 
thermochemical  matrix. 
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onooooooo  1000.  10.0  is.)  1.0  oct/air 

a 

1  HYDROGEN  1.008  Ifl. 

6  CARRON  12.011  ft. 

7  NITROGEN  14.00ft  -0.767 

ft  OXYGEN  16.000  -0.233 


1628000000000  OJANAF  01/61  C02 

-040540*5  165350*5  144550*2  2 1 03R6-3- 1 82302*7  798480*2  500.  3000.1  O.C02 

-940540*5  165350*5  156451 *2-381561 -4-602768*7  798480*2  3000.  5000.1  O.C02 

21  lftOOOOOOOOO  OJANAF  03/61  H20 

-577980*5  302010*5  112254*2  811397-3-260800*7  684210*2  500.  3000.1  0.H20 

-577980*5  302010*5  157278*2-191548-3-173599*8  684210*2  3000.  5000.1  0.H2O 

2700000000000  OJANAF  03/61  N2 

000000-0  221650*5  862699*1  116090-3-103715*7  637650*2  500.  3000.1  0.N2 

000000-0  221650*5  984175*1-116232-3-612728*7  637650*2  3000.  5000.1  0.N2 

1610000000000  OJANAF  03/61  CO 

-264170*5  223570*5  865040*1  117021-3-898211*6  653700*2  500.  3000.1  O.CO 

-264170*5  223570*5  115496*2-424139-3-131563*8  653700*2  3000.  5000.1  O.CO 

1718000000000  OJANAF  06/63  NO 

215800*5  227000*5  877623*1  899031-4-789656*6  688490*2  500.  3000.1  O.NO 

215800*5  227000*5  916260*1  657885-5-212519*7  688490*2  3000.  5000.1  O.NO 

1720000000000  OJANAF  06/63  N02 

801100*4  345800*5  137819*2  315611-4-133765*7  848890*2  500.  3000.1  0.NO2 

801100*4  345800*5  130154*2  172506-3  175340*7  848890*2  3000.  5000.1  0.NO2 

1738000000000  OJANAF  12/64  N03 

170000*5  498219*5  192301*2  191161-3-112129*7  998970*2  500.  3000.1  0.NO3 

170000*5  490219*5  198709*2  577142-7-172913*7  998970*2  3000.  SOOO.l  0.NO3 

2718000000000  OJANAF  12/60  N20 

195000*5  365450*5  144686*2  120489-3-153478*7  815900*2  500.  3000.1  0.N2O 

195000*5  365450*5  123036*2  459018-3  942302*7  015900*2  3000.  5000.1  O.N20 

2738000000000  OJANAF  12/64  N203 

197999*5  617799*5  232121*2  477580-3-124103*7  123217*3  500.  3000.1  0.N2O3 

197999*5  617799*5  240304*2  102872-5-293892*7  123217*3  3000.  5000.1  0.N2O3 

2748000000000  OJANAF  9/64  N204 

216999*4  785969*5  297734*2  594431-3-178033*7  134908*3  500.  3000.1  0.N2O4 

216999*4  785969*5  317015*2  175630-5-385969*7  134900*3  3000.  5000.1  O.N204 

2758000000000  OJANAF  12/64  N205 

216999*4  915759*5  353051*2  147295-3-163850*7  156719*3  500.  3000.1  0.N2O5 

216999*4  915759*5  357365*2  624545-5-171267*7  156719*3  3000.  5000.1  0.N205 

1  1  16  ‘0  00000000  OJANAF  03/61  CH 

142006*6  221300*5  026079*1  302211-3-100184*7  616120*2  500.  3000.1  O.CM 

142006*6  221300*5  707091*1  463281-3  552860*7  616120*2  3000.  SOOO.l  O.CM 

2116000000000  OJANAF  12/6?  CH2 

950000*5  329960*5  132894*2  413822-3-290273*7  684940*2  500.  3000.1  0.CH2 

950000*5  329960*5  140720*2  132150-3-239493*7  684940*2  3000.  5000.1  0.CH2 

3116000000000  OJANAF  12/62  CH3 

319400*5  434190*5  182763*2  401025-3-461203*7  786040*2  500.  3000.1  0.CH3 

319400*5  434190*5  204899*2-108028-3-1 14692*8  786040*2  3000.  5000.1  0.CH3 

4116000000000  OJANAF  03/61  CH4 

-178950*5  530790*5  230948*2  677896-3-755061*7  825970*2  500.  3000.1  0.CH4 

-178950*5  530790*5  236053*2  374323-3-368234*7  825970*2  3000.  5000.1  0.CM4 

1  126000000000  OOUFF  BAUER  6/61  C2H 

117395*6  349620*5  134210*2  469100-3-187509*7  781140*2  500.  3000.1  0.C2H 

117395*6  349620*5  148516*2  1 09402-3-S0386?*7  781140*2  3000.  SOOO.l  0.C2M 

2126000000000  OJANAF  0  3/61  C2H2 

541900*5  482570*5  109960*2  769044-3-409039*7  849690*2  500.  3000.1  0.C2H2 

541900*5  482570*5  203952*2  389062-3-645297*7  849690*2  3000.  5000.1  0.C2H2 

3126000000000  ODUFF  BAUER  6/61  C2H3 

659250*5  562430*5  230065*2  691322-3-499095*7  965150*2  500.  3000.1  0.C2H3 

659250*5  562430*5  239033*2  331771-3-335425*7  965150*2  3000.  5000.1  0.C2H3 

4126000000000  OJANAF  12/60  C2H4 

124960*5  676030*5  294887*2  526568-3-865313*7  101790*3  500.  3000.1  0.C2H4 

124960*5  676830*5  313872*2  616230-4-131058*8  101790*3  3000.  5000.1  0.C2M4 

6126000000000  OOUFF  BAUER  6/61  C2H6 

-203200*5  892710*5  404236*2-169721-2-193680*8  119447*3  500.  3000.1  0.C2H6 

-203200*5  892710*5  393119*2  764077-3-383891*7  119447*3  3000.  5000.1  0.C2M6 

1  136000000000  ODUFF  BAUER  6/61  C3H 

127703*6  409620*5  194464*2  508379-3-311933*7  928820*2  500.  3000.1  0.C3M 

127703*6  489620*5  199582*2  245687-3-632514*6  920820*2  3000.  5000.1  0.C3H 

2136000000000  ODUFF  BAUER  6/61  C3M2 

106522*6  635170*5  247444*2  992918-3-437596*7  104666*3  500.  3000.1  0.C3M2 

106522*6  635170*5  266623*2  364767-3-467730*7  104666*3  3000.  5000.1  0.C3H2 

3136000000000  ODUFF  BAUER  6/61  C3H3 

764850*5  709520**?  286929*2  766031-3-548098*7  113604*3  500.  3000.1  0.C3H3 

764050*5  709520*5  291672*2  456238-3-130601*7  113604*3  3000.  5000.1  0.C3H3 

4136000000000  OOUFF  BAUER  6/61  C3H4 

442940*5  822470*5  372219*2-528861-4-1 0 1994*0  120293*3  500.  3000.1  0.C3H4A 

442940*5  822479*5  344271*2  578173-3-200515*7  120293*3  3000.  5000.1  0.C3H4A 

51  360  00000000  ODUFF  BAUER  6/61  C3H5 

324310*5  946210*5  397325*2  958S55-3-904683*7  141759*3  500.  3000.1  0.C3H5 

324310*5  946210*5  39S073*?  733720-3-949766*6  141759*3  3000.  5000.1  0.C3H5 

1  146000000000  ODUFF  BAUER  6/61  C4H 

155196*6  645030*5  252172*2  574336-3-339547*7  113166*3  500.  3000.1  0.C4H 

155196*6  645030*5  249358*2  402513-3  377676*7  113166*3  3000.  5000.1  0.C4H 

2146000000000  OOUFF  BAUER  6/61  C4H2 

111715*6  763230*5  290097*2  735927-3-430646*7  119311*3  500.  3000.1  0.C4H2 

111715*6  763230*5  301707*2  456896-3  698124*6  119311*3  3000.  SOOO.l  0.C4H2 

3146000000000  ODUFF  BAUER  6/61  C4H3 

101975*6  851450*5  343217*2  855861-3-624935*7  130072*3  500.  3000.1  0.C4H3 

101975*6  851450*5  345194*2  564141-3-151557*6  130072*3  3000.  5000.1  0.C4H3 

4146000000000  ODUFF  BAUER  6/61  C4H4 

737050*5  9^1 ?40*5  398132*2  961427-3-792442*7  139440*3  500.  3000.1  0.C4H4A 

737050*5  97124**5  395899*2  720956-3  577880*6  139440*3  3000.  5000.1  0.C4H4A 

1  156000000000  OOUFF  BAUER  6/61  C5H 

185980*6  776470*5  299470*2  946114-3-430591*7  121859*3  500.  3000.1  0.C5H 

105980*6  776470*5  291522*2  690014-3  976229*7  121059*3  3000.  5000.1  0.C5H 


Sample  Problem  Number  1 
Listing  of  Input 
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?  1  s 

6  0  0 

o  o  o 

0000  ODUFF  BAUER  6/61 

C5H2 

lfc5?10*6 

933660*5 

169008*2 

965961-3-583890*7  137191*3 

500* 

3000.1 

0.C5H2 

165? 10*6 

933660*5 

370815*2 

618020-3  137008*7  137191*3 

3000* 

5000.1 

O.CSH? 

1  1  5 

6  0  0 

0  0  0 

0000  ODUFF  BAUER  6/61 

C5H3 

135414*6 

100517*6 

405177*2 

799736-3-675063*7  147599*3 

500* 

3000.1 

0.C5H3 

! 354 14 *6 

100517*6 

396951*2 

706570-3  334017*7  147599*3 

3000* 

5000.1 

0.C5H3 

1  1  6 

6  C  0 

0  0  0 

0000  ODUFF  BAUER  6/61 

C6H 

21 3164*6 

942310*5 

374077*2 

420402-3-474514*7  140238*3 

500* 

3000.1 

0.C6H 

?1 3164*6 

942310*5 

351470*2 

661760-3  962699*7  148238*3 

3000* 

5000.1 

0.C6H 

?  1  6 

6  0  0 

0  0  0 

0000  ODUFF  BAUER  6/61 

C6H2 

169663*6 

105052*6 

412157*2 

896814-3-569006*7  153063*3 

500. 

3000.1 

0.C6H2 

169663*6 

105052*6 

409176*2 

665591-3  323615*7  153063*3 

3000* 

5000.1 

0.C6H2 

3  1  6 

6  0  0 

ooo 

0000  ODUFF  BAUER  6/6! 

C6H3 

156461*6 

114639*6 

461927*2 

874250-3-738470*7  164599*3 

500* 

3000.1 

0.C6H3 

156461*6 

114639*6 

450580*2 

012263-3  450167*7  164599*3 

3000. 

5000.1 

0.C6H3 

6  1  6 

6  0  0 

0  0  0 

0000  ODUFF  BAUER  6/61 

C6H6 

197700*5 

147317*6 

626122*2 

117000-2-139277*0  172317*3 

500* 

3000.1 

0.C6H6 

197  700*5 

147317*6 

595572*2 

140760-2  737034*7  172317*3 

3000* 

5000.1 

0.C6H6 

1  1  7 

6  0  0 

0  0  0 

0000  ODUFF  BAUER  6/61 

C7H 

241752*6 

108509*6 

434324*2 

412320-3-559^94*7  160749*3 

500. 

3000.1 

0.C7M 

?4l 752*6 

108509*6 

402025*2 

023702-3  123670*0  160749*3 

3000* 

SOOO.l 

0.C7H 

2  1  7 

6  0  0 

0  0  0 

0000  OOUFF  BAUER  6/61 

C7M2 

220830 '6 

i 19071*6 

479372*2 

119639-2-020430*7  166034*3 

500* 

3000.1 

0.C7H2 

220030*6 

119871*6 

479300*2 

025569-3  171374*7  166034*3 

3000* 

5000.1 

0.C7H2 

1  1  8 

6  0  0 

0  0  0 

0000  ODUFF  BAUER  6/61 

C8H 

?88675*6 

122402*6 

489007*2 

527442-3-644448*7  174023*3 

500* 

3000.1 

O.CSH 

200075*6 

122402*6 

454317*2 

953953-3  140527*0  174023*3 

3000* 

5000.1 

O.CSH 

2  1  6 

6  0  0 

0  0  0 

0000  ODUFF  BAUER  6/61 

CSH2 

226156*6 

133849*6 

520472*2 

975001-3-710651*7  100209*3 

500* 

3000.1 

O.CSH? 

226156*6 

133049*6 

511890*2 

960001-3  021940*7  100289*3 

3000* 

5000.1 

O.C0H2 

1  1  9 

6  0  0 

0  0  0 

0000  ODUFF  BAUER  6/61 

C9H 

293122*6 

137300*6 

546536*2 

609606-3-683510*7  189903*3 

500* 

3000.1 

0.C9H 

293122*6 

137300*6 

507642*2 

106604-2  150457*0  109903*3 

3000* 

5000.1 

0.C9H 

2  1  9 

6  0  0 

0  0  0 

0  0  •  0  ODUFF  BAUER  6/6! 

C9H2 

272220*6 

140147*6 

*92585*2 

135051-2-900554*7  194630*3 

500* 

3000.1 

0.C9H2 

272?20*6 

140147*6 

505236*2 

106325-2  470277*7  194630*3 

3000* 

5000.1 

0.C9H2 

1  1  10 

6  0  0 

0  0  0 

0000  ODUFF  BAUER  6/61 

C10H 

346750*6 

151795*6 

605449*2 

627977-3-760351*7  208392*3 

500. 

3000.1 

0.C10H 

346750*6 

151795*6 

562113*2 

115764-2  170972*0  200392*3 

3000* 

5000.1 

O.CIOH 

2  1  10 

6  0  0 

0  0  0 

0000  OOUFF  BAUER  6/61 

C10H2 

264C31 *6 

163004*6 

643250*2 

10<m5-2-83184S»7  214369*3 

500. 

3000.1 

0.C10H2 

284031*6 

163004*6 

620020*2 

11S83Z-2  109441*6  214369*3 

3000* 

5000.1 

0.C10H2 

1  1  1 

6  l  0 

0  0  0 

0000  0J4N4F  03/61 

CHO 

■2900U0*4 

323670*5 

120033*2 

300636-3-201721*7  789830*2 

500* 

3000.1 

O.CHO 

■290000*4 

123670*5 

103020*2 

633312-3  I216I5**  789830*2 

3000* 

5000.1 

O.CHO 

2  l  1 

6  1  8 

0  0  0 

0000  0J6N4F  03/61 

CH20 

■277000*5 

437910*5 

184022*2 

379921-3-431962*7  848880*2 

500. 

3000.1 

0.CH2O 

•277000*5 

437910*5 

10*620*2 

161963-3-372770*7  848880*2 

3000. 

5000.1 

0.CH2O 

4  1  2 

6  1  8 

0  0  0 

0000  0J4N4F  12/60 

C2H40 

•121900*5 

013210*5 

353747*2 

546144-3-964936*7  117211*3 

500* 

3000.1 

O.C2H40 

•121900*5 

013210*5 

373112*2 

690621-4-141966*8  117211*3 

3000* 

5000.1 

O.C2H40 

1  1  1 

6  1  7 

0  0  0 

0000  OJANAF  03/61 

CHN 

312000*5 

155930*5 

137021*2 

552243-3-228955*7  758620*2 

500* 

3000.1 

O.CHN 

312000*5 

355930*5 

170895*2-295052-3-103671*8  758620*2 

3000* 

5000.1 

O.CHN 

1  1  1 

6  1  7 

1  0  0 

0000  OJANAF  12/60 

CHNO 

•279000*5 

461930*5 

179701*2 

469024-3-200904*7  920990*2 

500. 

3000.1 

O.CHNO 

■279000*5 

461930*5 

210167*2-229915-3-105560*0  928990*2 

3000. 

5000.1 

O.CHNO 

1  1  0 

0  0  0 

0  0  0 

0000  OJANAF  12/60 

H 

521020*5 

134230*5 

400223*1 

555469-4  173095*6  306620*2 

500. 

3000.1 

O.H 

521020*5 

134230*5 

300752*1 

3150*5-3  929744*7  380620*2 

3000. 

5000.1 

O.H 

1  7  0 

0  0  0 

0  0  0 

0000  OJANAF  01/61 

N 

112965*6 

134370*5 

486944* 1 

303516-4  950460*5  480900*2 

500. 

3000.1 

O.N 

11 2965*6 

134370*5 

420957*1 

240844-3-417273*6  400900*2 

3000. 

5000.1 

O.N 

1  6  0 

0  0  0 

0  0  0 

0000  OJANAF  06/62 

0 

595593*5 

135220*5 

497228*1 

300760-5  154749*5  500960*2 

500. 

3000.1 

0.0 

595590*5 

135220*5 

657489*1-224260-3-891782*7  500960*2 

3000. 

5000,1 

0.0 

2  8  0 

0  0  0 

0  0  0 

0000  OJANAF  03/61 

02 

000000-0 

234460*5 

804370*1 

510072-3-152710*6  679730*2 

500. 

3000.1 

0.02 

000000-0 

234460*5 

101071*2 

290991-4-703079*7  679730*2 

3000. 

5000.1 

0.02 

3  H  0 

0  0  0 

0  0  0 

0000  OJANAF  6/61 

03 

341000*5 

160230*5 

135311*2 

230739-3-606511*6  860919*2 

500. 

3000.1 

0.03 

341000*5 

160230*5 

139121*2 

125613-3-900654*6  860919*2 

3000. 

5000.1 

0.03 

1  1  1 

7  0  0 

0  0  0 

0000  OJANAF  12/60 

HN 

792r00*s 

217660*5 

8211 31*1 

2815550-126096*7  609290*2 

500. 

3000.1 

O.MN 

792000*5 

217660*5 

765457*1 

359053-3  180770*7  609290*2 

3000. 

SOOO.l 

0  .HN 

1  1  1 

7  1  8 

0  0  0 

0000  OJANAF  12/60 

HNO 

236000*5 

327550*5 

133820*2 

122061-3-232403*7  703050*2 

500. 

3000.1 

O.HNO 

236000*S 

327550*5 

141523*2-670424-4-4691 35*7  783850*2 

3000. 

5000.1 

O.HNO 

1  1  1 

7  2  8 

0  0  0 

0000  OJANAF  06/63 

HN02 

■183400*5 

468800*5 

190270*2 

189744-3-300752*7  962310*2 

500. 

3000.1 

O.HNO 2 

•103400*5 

460800*5 

210284*2- 

■229017-3-970618*7  962310*2 

3000. 

5000,1 

0.HNO2 

1  1  1 

7  18 

0  0  0 

0000  OJANAF  06/63 

HNO  3 

■321000*5 

609120*5 

253548*2 

600154-4-450532*7  110051*3 

500. 

3000.1 

0.HN03 

•321000*5 

609120*5 

266283*2-145705-3-102656*0  110851*3 

3000. 

5000.1 

0.HNO3 

1  1  1 

8  0  0 

0  0  0 

0000  OJANAF  12/60 

HO 

933000*4 

214040*5 

773193*1 

394306**3-973561  *6  613820*2 

500. 

3000.1 

O.HO 

933000*4 

214040*5 

965144* 1 -443520-4-6861 15*7  613020*2 

3000. 

5000.1 

O.HO 

1  1  ? 

8  0  0 

0  0  0 

0000  OJANAF  03/64 

111 

M02 

500000*4 

325340*5 

127161*2 

292770-3-172279*7  799790*2 

1000. 

3000.1 

O.HO? 

500000*4 

325340*5 

138020*2 

147000-4-399551*7  799790*2 

3000. 

5000.1 

O.HO? 

2  l  0 

0  0  0 

0  0  0 

0000  OJANAF  03/61 

H2 

000000-0 

212100*5 

711963*1 

621 950- 3-71 2t» 94*6  404650*2 

500. 

3000.1 

0.H2 

000000-0 

212100*5 

681794*1 

509054-3  265106*7  404650*2 

3000. 

5000.1 

0.H2 

2  1  1 

7  0  0 

0  0  0 

0000  OJANAF  12/65 

H2N 

400999*5 

105799*5 

975862*1 

114401-2-510970*6  701500*2 

500. 

3000.1 

0.H2N 

400999*5 

105799*5 

117419*2 

229493-4-610024*7  701500*2 

3000. 

5000.1 

0.H2N 

Sample  Problem  Number  1 
Listing  of  Input,  Continued 
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I 


31  17000000000  OJANAF  09/65 

-109700*5  416309*5  125186*2  219670-2-978531*6  769080*2  500* 

-109700*5  416309*5  166320*2  923148-3*361337*7  769080*2  3000. 
4127000000000  OJANAF  12/65 
227899*5  685359*5  229856*2  241933-2-195931*7  108099*3  500. 

227899*5  685359*5  314319*2  487945-4-1 39705*8  108099*3  3000. 
3628000000000  OJANAF  12/60 
-830000*4  652120*5  263775*2  103904-3-350368*7  112723*3  500. 

-830000*4  652120*5  263543*2  88220 l -4-287091 *7  112723*3  3000. 
1617000000000  OJANAF  12/62 
109000*6  232490*5  655906*1  115326-2  479517*6  669760*2  500. 

109000*6  232490*5  988013*1  31 3855-3-649453* f  669760*2  3000. 
2627000000000  OJANAF  3/61 
738699*5  511070*5  188740*2  559856-3-896873*6  985479*2  500. 

738699*5  511070*5  208204*2  630229-5-346865*7  985479*2  3000. 
4627000000000  OJANAF  3/61 
127500*6  801239*3  29V490*2  856052-3-152255*7  133063*3  500. 

127500*6  801239*5  327212*2  952999-5-541654*7  133063*3  3000. 

1600000000000  OJANAF  0  3/61 
170886*6  135500*5  444433*1  228125-3  409830*6  492870*2  500. 

170886*6  135500*5  4122t2*l  261908-3  262886*7  492870*2  3000. 
2600000000000  OJANAF  9/61 
198999*6  246990*5  776612*1  696081-3  185649*6  685519*2  500. 

198999*6  246990*5  104162*2  566841-4-640205*7  685519*2  3000. 

3600000000000  OJANAF  12/60 
189670*6  366220*5  146441*2  622S36-4-168227*?  798410*2  500. 

189670*6  366220*5  144782*2  792232-4-646877*6  798410*2  3000. 
4600000000000  OJANAF  12/60 
242321*6  511239*5  205903*2  623436-4-257703*7  986760*2  500. 

242321*6  511230*5  210714*2-434895-4-404939*7  986760*2  3000. 

5600000000000  OJANAF  12/60 
242374*6  656230*5  264706*2  806528-4-337125*7  111641*3  500. 

242374*6  656230*5  271156*2-580271-4-543183*7  111641*3  3000. 

6600000000000  OOUFF  BAUER  6/61 
291342*6  772700*5  303965*2  607305-3-475660*7  128958*3  500. 

291342*6  772700*5  308488*2  643430-4  583350*7  128958*3  3000. 
7600000000000  OOUFF  BAUER  6/61 
292116*6  909770*5  359795*2  669838-3-579638*7  160555*3  500. 
292116*6  909770*5  355785*2  210202-3  102222*8  140555*3  3000. 

6600000000000  OOUFF  BAUER  6/61 
344891*6  106684*6  415594*2  733634-3-683462*7  159275*3  500. 

344891*6  ! 04684*6  402894*2  359448-3  146984*8  159275*3  3000. 

9600  000000000  OOUFF  BAUER  6/61 

340665*6  118391*6  471410*2  796582-3-787328*7  170872*3  500. 

340665*6  118391*6  450207*2  505833-3  190813*8  170872*3  3000. 
10  600000000000  OOUFF  BAUER  6/61 

400439*6  132098*6  527202*2  860613-3-891076*7  189592*3  500. 

400439*6  132098*6  497419*2  651908-3  235109*8  189592*3  3000. 

1600000000000  OJANAF  03/6T 
000000-0  144120*5  586075*1  953976-4-766621*6  121290*2  500. 

000000-0  144120*5  485114*1  291605-3  307202*7  121290*2  3000. 


H3N 

3004.1 

0.H3N 

5000.1 

0.H3N 

H4N2 

3000.1 

0.H4N2 

5000.1 

0.H4N2 

C302 

3900.1 

0.C3O2 

5000.1 

0.C3O2 

CN 

3000.1 

O.CN 

5000.1 

O.CN 

C2N2 

3000.1 

0.C2N2 

S000.1 

0.C2N2 

C4N2 

3000.1 

0.C4N2 

5000.1 

0.C4N2 

C 

3000.1 

o.c 

5000.1 

o.c 

C2 

3000.1 

O.C? 

5000.1 

O.C? 

C3 

3090.1 

0.C3 

5000.1 

0.C3 

C4 

3000.1 

0.C4 

5000.1 

0.C4 

C5 

3000.1 

0.C5 

5000.1 

0  .C5 

C6 

3000.1 

0.C6 

5000.1 

0 .06 

C7 

3000.1 

0.C7 

5000.1 

0.C7 

08 

3000.1 

0.C8 

5000.1 

0.C8 

C9 

3000.1 

0.C9 

5000.1 

0.C9 

CIO 

3000.1 

0.C10 

5000.1 

0.C10 

C« 

3000.2 

100. C* 

5000.2 

100. c* 

2500. 

2000. 

500. 

1000. 

2500. 

2000. 

500. 


OCT/ AIR 
OCT/ AIR 
OCT/AIR 

16.5  1.0  OCT/AIR 

OCT/AIR 

OCT/AIR 

OCT/AIR 


Sample  Problem  Number  1 
Lifting  of  Input,  Concluded 
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AEROThCRM 

EQUILIBRIUM 

SURFACE  THEWMOCHf histry 

JOd  OCT/AIR 

Oil QOOOOO 

relative  ELEMENTAL  COMPOSITIONS,  atomic  a  I  s/unit 

HASS 

AT, no. 

ELEMENT 

ATOMIC 

COMPONENT 

1  COMPONENT  2 

COMPONENT  3 

1 

HYOHOGf N 

1.00000 

.0000000 

.1575741 

-.0000000 

6 

CARBON 

12,01100 

.0000000  ' 

.0700329 

-.0000000 

7 

NITROGEN 

14.00000 

•  054  7544 

.0000000 

-.0000000 

0 

OXYGEN 

16.00000 

.0145625 

.0000000 

-.0000000 

element 

HYDROGEN 

CARBON 

NITROGEN 

OXYGEN 

BASE  SP 

H20 

C02 

N2 

CO 

AEROThERM  EQUILIBRIUM  SURFACE  THERMOCHEMISTRY  SOLUTION  OCT/AIR 


HcRIVATlVE  PROPERTY  OUTPUT 

CP-FROZEN  CP-EQUIL  DLNM  ULNT  OLNH/OLNP  GAMMA 
.35365-00  .81113-00  -.20037-00  .12870-01  .11543*01 

PROPERTY  ROUTINE  OUTPUT  IN  LB-MASS. F T .SEC.0TU. AND  DEG-R 

MU1  MU2  M0L.8T  HTIL  CPTIL  HTIL* 

.10590*01  .25361*0?  .27406*02  .83312*03  .37030-00  .77023*03 

ELEMENTAL  K  AND  Z  MASS  FRACTIONS  BY  ATOMIC  NUMBER  ....  (GAMEX  »  ,6671 

16  7  8 

.96655-02  .52246-01  .71936-00  .21053-00 

.12351-01  .*8658-01  .72005-00  .21094-00 


SOLUTION  TIMES 

ITERATIONS  *  46 


MACMlNt  TIME  *  3.631  SEC. 


CLOSED  SYSTEM  EQUILIBRIUM  SOLUTION  OUTPUT 

RELATIVE  MASSES  OF  COMPONENTS  1.2.  AND  3.  .15100*02  .10000*01  -.00000 

THERMODYNAMIC  STATE 


PRESSURE 

(ATMOSPHERES) 

a 

10.0000 

TEMPERA  JURE  (OEG  Kl 

■  3000.000 

(LBF/IN2-ABS) 

a 

146.9600 

(OEG  R) 

*  5400.000 

OENSITY 

(GM/CM3) 

s 

.1113-02 

ENTHALPY 

(CAL/GM) 

•  394.846 

(LBM/FT3) 

* 

.6949-01 

(BTU/LBM) 

»  710.726 

MOLECULAR  HEIGHT 

3 

27.4063 

ENTROPY 

(0TU/L8M-R) 

*  .2295*01 

CHEMICAL  STATE 

MOLE  FR. 

*  MOLECULES  / 

TOTAL  GAS  PHASE  MOLECULES 

SPECIES 

MOLE  FR. 

SPECIES  MOLE  FR. 

SPECIES 

MOLE  FR. 

SPECIES 

MOLE  FR. 

SPECIES 

MOLE  FR* 

C02 

.60344-01 

H20 

.11107*00 

N2 

•69661-00 

CO 

*50869-01 

NO 

.14174*01 

N02 

.70660-05 

NO  3 

.27093-10 

N20 

•21844-05 

N203 

•51878-12 

N204 

•25019-16 

N205 

•19441-19 

CH 

.36906-12 

CH2 

.51277-13 

CH3 

•57583-12 

CH4 

•36782-13 

C2H 

.12183-15 

C2H2 

.16826-15 

C2H3 

.20809-19 

C2H4 

.34354-20 

C2H6 

•16763-25 

C3H 

.39459-22 

C3M2 

.45716-24 

C3H3 

.18447-25 

C3H4 

.18026-27 

C3HS 

•82435-28 

C4H 

.80602-29 

C4H2 

.40108-30 

C4H3 

.10658-32 

C4H4 

•16523-34 

C5H 

•50215-38 

CSH2 

.00000 

CSH3 

•00000 

C6H 

•00000 

C6H2 

•00000 

C6H3 

•00000 

C6H6 

.00000 

C7H 

•OOOOO 

C7H2 

•OOOOO 

C8H 

•00000 

C8H2 

•00000 

C9M 

.00000 

C9H2 

•00000 

C10H 

•OOOOO 

C10H2 

•OOOOO 

CHO 

•17564-05 

CH20 

.32759-08 

C2M40 

.60298-23 

CHN 

•10623-07 

CHNO 

•22941-07 

H 

•53175*02 

N 

.36473-05 

0 

.49424-02 

02 

.193*3-01 

03 

•82171*06 

HN 

•16865*05 

HNO 

.22230-05 

HN02 

.22778-06 

HN03 

.39993-10 

HO 

•17850-01 

H02 

•77561*05 

H2 

.11453-01 

H2N 

.28597-06 

H3N 

.71272-07 

H4N2 

•32820-15 

C302 

•61663*18 

CN 

.20774-09 

C2N2 

.11532-15 

C4N2 

.49147-30 

C 

•24470-11 

C2 

•62000-16 

C3 

.13318-23 

C4 

.25244-32 

C5 

.16951-38 

C6 

•00000 

C? 

•OOOOO 

C8 

.00000 

C9 

•OOOOO 

CIO 

•00000 

C* 

•00000 

Sample  Problem  Number  1 
Sample  Output 
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AEROTHERM  EQUILIBRIUM  SURFACE  THERMOCHEMISTRY1  JOB  OCT/AJW 


DERIVATIVE  PROPERTY  OUTPUT 

CP-FROZEN  CP-EQU1L  DLNM/OLNT  OLNM/DLNP  GAMMA 
.35035-00  .7*342-00  -.27215-00  .12526-01  . 11553*01 


PROPERTY  ROUTINE  OUTPUT  IN  LB-MmSS*FT.SEC»BTU,AND  DEG-R 

MU1  MU2  MOL. HI  HT1L  CRT IL 

.10637*01  .25473*02  *275  5*02  .86527*03  .36474-00 


HTIL* 

.81653*03 


elemental  K  ANO  Z  MASS  fractions  by  atomic  number 

16  7  6 

.90763-02  .48067-01  .72317-00  .21969-00 

.11243-01  .44563-01  .72461-00  .21958-00 


1GAMEX  =  .667) 


SOLUTION  TIMES 

ITERATIONS  *  12 


MACHINE  TIME  »  .611  SEC. 


CLOSED  SYSTEM  EQUILIBRIUM  SOLUTION  OUTPUT 


RELATIVE  MASSES  OF  COMPONENTS 

1*2.  ANO  3. 

•16500*02  .10000*01 

-.00000 

THERMOOYNAMIC  STATE 
PRESSURE  (ATMOSPHERES) 

E 

10.0000 

TEMPERATURE  (DEG  K) 

X 

3000.000 

(L8F/IN2-ABS) 

= 

146.9600 

(DEG  R) 

X 

5400.000 

OENSITY  (GM/CM3) 

- 

.1120-02 

ENTHALPY  (CAL/GM) 

X 

419,191 

(LBM/FT3) 

= 

.6992-01 

(BTU/LBM) 

9 

754,543 

MOLECULAR  HEIGHT 

S 

27.5750 

ENTROPY  (BTU/LBM-R) 

X 

.2282*01 

CHEMICAL  STATE 

MOLE  PR.  »  MOLECULES  /  TOTAL  GAS  PHASE  MOLECULES 


SPECIES 

MOLE  FR. 

SPECIES 

MOLE  FR. 

SPECIES 

MOLE  FR. 

SPECIES 

MOLE  FR. 

SPECIES 

MOLE  FR. 

C02 

.67964-01 

H20 

.10344*00 

N2 

.70328-00 

CO 

•42387-01 

NO 

.16995-01 

N02 

.10113-04 

N03 

.46272-10 

N20 

.26319-05 

N203 

•89026-12 

N204 

•51240-16 

N205 

.47517-19 

CH 

.22763-12 

CH2 

.27939-13 

CH3 

•27715-12 

CH4 

•15639-13 

C2H 

.52467-16 

C2H2 

.64009-16 

C2H3 

.69929-20 

C2M4 

•10199-20 

C2M6 

.38833-26 

C3H 

.11864-22 

C3H2 

.12143-24 

C3H3 

.43283-26 

C3H4 

•37363-28 

C3HS 

•15094-20 

C4H 

.18617-29 

C4H2 

.74381-31 

C4H3 

.17460-33 

C4M4 

.26806-35 

C5H 

•00000 

C5H2 

•00000 

C5M3 

.00000 

C6H 

.00000 

C6M2 

•00000 

C6H3 

•00000 

C6H6 

.00000 

C7H 

•00000 

C7H2 

.00000 

C8H 

•00000 

C8M2 

•00000 

C9H 

.00000 

C9H2 

.0U000 

C10H 

•00000 

C10H2 

•00000 

CHO 

.12929-05 

CH20 

.21301-08 

C2H40 

.21363-23 

CHN 

.67071-08 

CHNO 

•16968-07 

H 

.46974-02 

N 

.36647-05 

0 

.58984-02 

02 

.27550-01 

03 

.13967-07 

HN 

•14970-05 

HNO 

.23548-05 

hno2 

.28796-06 

HN03 

.60339-10 

HO 

•18618-01 

H02 

•97611-05 

H2 

.89376-02 

H2N 

.22423-06 

H3N 

.49366-07 

H4N2 

.20177-15 

C302 

.29893-10 

CN 

.14574-09 

C2N2 

.56756-16 

C4N2 

•11792-30 

C 

•17085-11 

C2 

.30224-10 

C3 

.45329-24 

C4 

.59992-33 

C5 

•00000 

C6 

•00000 

C7 

•OOOOO 

C8 

.00000 

C9 

•00000 

CIO 

•00000 

C* 

•ooooo 

Sample  Problem  Number  1 
Sample  Output,  Concluded 
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Ptl  1  1  -777.4  6R.  .70*  .0*0  .10  AL  PROPLNT 

COMPOSITION  or  THE  BOUNDARY  LAYER  EDGE  FLOW  IN  A  SOLID  PROPELLANT  ROCKET  MOTOR 
PROPELLANT  IS  70.4  PCT  AMMONIUM  PERCHLORATE.  0.6  PCT  BUTADIENE  UNDER. 

ANO  ?0.0  PCT  ALUMINUM  POMDER.  PRESSURE  IS  66.  ATM.  ENTHALPY  IS  -TTT.  CAL/OM 


1  HYDROGEN 

1*008 

4.0 

6  CARBON 

I2.nn 

7  NITROGEN 

14.00a 

1.0 

0  OXYGEN 

16*000 

4.0 

13  ALUMINUM 

26.970 

17  CHLORINE 

35.457 

1.0 

1  1  o 

0  0  0 

0  0  0 

0000  OJANAF  12/60 

H 

521020*5 

134230*5 

480223*1 

555469-4  173095*6  388620*2 

500* 

3000*1 

O.H 

521020*5 

134230*5 

308752*1 

315025-3  929744*7  308620*2 

3000. 

5000.1 

O.H 

1  6  1 

8  0  0 

0  0  0 

0000  OJANAF  03/61 

CO 

-264170*5 

223570*5 

865040*1 

117021-3-898211*6  653700*2 

500* 

3000.1 

O.CO 

-264170*5 

223570*5 

115496*2-424139-3-131563*8  653700*2 

3000* 

5000.1 

O.CO 

1  7  0 

0  0  0 

0  0  0 

0000  OJANAF  03/61 

N 

1 1 2965*6 

134370*5 

486944*1 

3R3516-4  958460*5  480900*2 

500. 

3000.1 

O.N 

112965*6 

134370*5 

428957*1 

240844-3-417273*6  480900*2 

3000* 

5000*1 

O.N 

1  8  0 

0  0  0 

0  0  0 

0000  OJANAF  06/62 

0 

595590*5 

1 3S2?0*5 

497228*1 

380760-5  154749*5  500960*2 

500. 

3000.1 

0.0 

595590*5 

135220*5 

6S7489* 1-224268-3-691782*7  500960*2 

3000* 

5000.1 

0.0 

1  13  0 

0  0  0 

0  0  0 

0000  OJANAF  03/61 

A1 

780000*5 

134610*5 

492905*1 

113634-4  706967*5  508440*2 

500* 

3000.1 

O.AL 

780000*5 

134610*5 

395674*1 

204094-3  361777*7  500440*2 

3000* 

5000.1 

O.AL 

1001  01 

17 

JANAF  03/61 

CLH 

-219700*5 

218040*5 

839693*1 

214028-3-125475*7  623640*2 

500* 

3000.1 

O.CLM 

-219700*5 

218640*5 

837177*1 

ln326 1-3-446002*6  623640*2 

3000* 

5000.1 

O.CLH 

1  6  1 

13  0  0 

0  0  0 

0000  OJANAF  12/60 

ALC 

209958*6 

240440*5 

896909*1-757687-5-752532*5  746650*2 

500. 

3000.1 

O.ALC 

209950*6 

240440*5 

961609*1-128042-3-264567*7  746650*2 

3000* 

5000.1 

O.ALC 

1  13  1 

17  0  0 

0  0  0 

0  0  0  0  OJANAF  06/61 

ALCL 

-116200*5 

244690*5 

895305*1 

115091-3-929251*5  750380*2 

500* 

3000.1 

O.ALCL 

-116200*5 

244690*5 

951 009* l 

835109-5-223156*7  750380*2 

3000. 

5000.1 

O.ALCL 

1  13  ? 

17  0  0 

0  0  0 

0000  OJANAF  06/61 

ALCL  2 

-780000*5 

371170*5 

138896*2 

586012-5-137329*6  106813*3 

500* 

3000.1 

O.AL CL 2 

-780000*5 

371170*5 

136493*2 

486738-4  869684*6  100813*3 

3000. 

5000.1 

0.ALCL2 

1  8  2 

13  0  0 

0  0  0 

0000  OJANAF  09/61 

AL20 

-314400*5 

359700*5 

138091*2 

405413-5-663265*6  912008*2 

500. 

3000.1 

0.AL2O 

-314400*5 

3S9700*5 

1446S6*2~1 08570*3“278920*7  912000.2 

3000. 

5000.1 

O.AL 20 

1  6  0 

0  0  0 

0  0  0 

0  0  0  0  OJANAF  03/61 

C 

170886*6 

135500*5 

444433*1 

PPM  25-3  609830*6  4928T0.2 

500. 

3000.1 

o.c 

170806*6 

135500*5 

412212*1 

261908-3  262886. T  492870*2 

3000* 

5000.1 

O.C 

4  1  1 

6  0  0 

0  0  0 

0000  OJANAF  03/61 

CH4 

-178950*5 

530790*5 

230948*2 

677896-3-755061*7  825970*2 

500* 

3000.1 

0.CH4 

-170950*5 

530790*5 

236053*2 

374323-3-368234*7  825970*2 

3000. 

5000.1 

0.CH4 

1  6  2 

8  0  0 

0  0  0 

0000  OJANAF  03/61 

C02 

-940540*5 

365350*5 

144559*2 

210386-3-182392*7  798480*2 

500. 

3000.1 

0.CO2 

-940540*5 

3653S0*5 

156451 *2-381561-4-602768*7  796400*2 

3000* 

5000.1 

0.C02 

2  1  2 

6  0  0 

0  8  0 

0  0  o  o  OJANAF  03/61 

C2H2 

541900*5 

482570*5 

189960*2 

769044-3-409039*7  849690*2 

500. 

3000.1 

0.C2H2 

541900*5 

482570*5  203052*2 

309002-3-64^297*7  849690*2 

3000* 

5000.1 

0.C2H2 

3  6  0 

0  0  0 

0  0  0 

0000  OJANAF  12/60 

C3 

189670*6 

366220*5 

146441*2 

622536-4-168227*7  790410*2 

500. 

3000.1 

0.C3 

109670*6 

366220*5 

144782*2 

792232-4-646077*6  790410*2 

3000. 

5000.1 

0.C3 

1  17  0 

0  0  0 

0  0  0 

0000  OJANAF  03/61 

a 

289220*5 

uo?so*$ 

508338*1-272037-4  290013*6  516230*2 

500. 

3000.) 

O.CL 

289220*5 

140280*5 

539251*1-770696-4-1 14573*7  516230*2 

3000* 

5000.1 

O.CL 

2  17  0 

0  0  0 

0  0  0 

0000  OJANAF  03/61 

CL2 

000000-0 

244290*5 

852234*1 

273227-3  207718*6  737600*2 

500* 

3000.1 

0.CL2 

000000-0 

244290*5 

100329*2-141 907-4-554731 ♦ 7  737600*2 

3000* 

5000.1 

O.CL  2 

1  1  1 

BOO 

0  0  0 

0000  OJANAF  12/60 

HO 

933000*4 

214040*5 

773193*1 

394306-3-973561*6  613820*2 

soo. 

3000.1 

0  .HO 

933000*4 

214040*5 

965144*1-443528-4-686115*7  613820*2 

3000* 

5000.1 

O.HO 

2  1  0 

0  0  0 

ooo 

0000  OJANAF  03/61 

H2 

000000-0 

212100*5 

711963*1 

621990-3— 712694*6  484650*2 

500* 

3000.1 

0.H2 

oooooo-o 

212100*5 

681794*1 

589054-3  265106*7  484650*2 

3000. 

5000.1 

0.H2 

2  1  1 

8  0  0 

0  0  0 

0000  OJANAF  03/61 

H20 

-577980*5 

302010*5 

112254*2 

811397-3-260800*7  684210*2 

500. 

3000.1 

O.H20 

-577900*5 

302010*5 

157270*2-191548-3-173599*8  604210*2 

3000. 

5000.1 

O.H20 

1  7  1 

0  0  0 

0  0  0 

0000  OJANAF  06/63 

NO 

215800*5 

227000*5 

077623*1 

8990 3 1 -4-709656*6  608490*2 

500. 

3000.1 

O.NO 

215800*5 

227000*5 

916260*1 

657805-5-212519*7  688490*2 

3000. 

5000.1 

0  .NO 

2  7  0 

0  0  0 

0  0  0 

0000  OJANAF  03/61 

N2 

000000-0 

221650*5 

862699*1 

116090-3-103715*7  637650*2 

500. 

3000.1 

O.N? 

000000-0 

221650*5 

984 175*1-11 6232-3-612728*7  637650*2 

3000. 

5000.1 

0.N2 

2  8  0 

0  0  0 

0  0  0 

0000  OJANAF  03/6! 

02 

000000-0 

234460*5 

804370*1 

510872r3-152718*6  679730*2 

500. 

3000.1 

0.02 

000000-0 

234460*9 

103071*2 

29099144-783079*7  679730*2 

3000. 

5000.1 

0*02 

3  8  2 

13  0  0 

0  0  0 

0000  OJANAF  09/61 

AL203* 

-400400*6 

811817*5 

319699*2 

120451-4-217514*7  772577*2 

500. 

2316.2 

100.AL2O3* 

-366770*6 

759318*5 

351 153*2-192777-4-555942*6  893107*2 

2310. 

5000.3 

100. At 203* 

1  6  0 

0  0  0 

0  0  0 

0000  OJANAF  03/61 

C* 

000000-0 

144120*5 

586075*1 

953976-4-766621*6  121290*2 

500. 

3000.2 

loo.c* 

000000-0 

144120*5 

405134*1 

291605-3  307202*7  121290*2 

3000* 

5000,2 

100. C* 

?  -1360.1  ?.00 

R 


Sample  Problem  Number  2 
Listing  of  Input 
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AEROTHERM  equilibrium  surface  thermochemistry  JOB  AL  PROPLNT  2U*1«1«* 
COMPOSITION  OF  THE  BOUNOARY  LAYER  EDGE  FLO*  IN  A  SOLID  PROPELLANT  ROCKET  MOTOR 
PROPELLANT  IS  70.4  PCT  AMMONIUM  PERCHLORATE ♦  9.t>  PCT  BUTADIENE  BINDER* 

AND  20.0  PCT  ALUMINUM  POWDER.  PRESSURE  IS  68.  ATM*  ENTHALPY  IS  -777.  CAL/GM 

RELATIVE  ELEMENTAL  COMPOSITIONS.  ATOMIC  WTS/UNIT  MASS 


AT. NO. 

element 

ATOMIC  WT 

COMPONENT  1 

l  COMPONENT 

2  COMPONENT 

3 

1 

HYDROGEN 

1.00800 

.0340434 

.1109221 

•0000000 

6 

CARBON 

12.01100 

•0000000 

.0739481 

.0000000 

7 

NITROGEN 

14.00800 

.0085109 

.0000000 

.0000000 

e 

OXYGEN 

16.00000 

.0340434 

.0000000 

•0000000 

13 

ALUMINUM 

26.97000 

•0000000 

.0000000 

.0370782 

17 

CHLORINE 

35.45700 

.0085109 

•00f 0000 

.0000000 

ELEMENT 

HYOROGEN 

CARBON 

NITROGEN 

OXYGEN 

ALUMINUM 

CHLORINE 

BASE  SP 

H 

CO 

N 

0 

AL 

CLH 

ITS 

i  TEMP 

PRES*MWT 

EQUIL  ER 

masbal  e« 

SCALE 

1  MASBAL 

2  MASBAL 

3  MASBAL 

4  MASBAL 

5  MASBAL 

6  MASBAL 

1 

3000.00 

.136*04 

.608*02 

.259*01 

.000 

•212*02 

-.199*02 

-.368*01 

.229*02 

-.765*01 

-.125*02 

2 

3297.98 

•130*04 

•506*02 

.276*01 

.169-00 

.172*02 

-.161*02 

-.297*01 

•186*02 

-.620*01 

-.102*02 

3 

3493.82 

.128*04 

.362*02 

.340*01 

.174-00 

•142*02 

-.133*02 

-.246*01 

•154*02 

-.512*01 

-.840*01 

4 

3582.42 

.129*04 

.282*02 

.518*01 

.136-00 

.123*02 

-.115*02 

-.212*01 

•133*02 

-.442*01 

-.726*01 

S 

3603.64 

.133*04 

.203*02 

.449*01 

.702-01 

.115*02 

-.107*02 

-.197*01 

•124*02 

-.411*01 

-.675*01 

6 

3631.40 

•140*04 

.185*02 

.316*01 

.850-01 

.105*02 

-.979*01 

-.179*01 

•114*02 

-.374*01 

-.616*01 

7 

3667.30 

.151*04 

.162*02 

.142*01 

.126-00 

.933*01 

-.852*01 

-.154*01 

•*00*02 

-.324*01 

-.536*01 

a 

3703.21 

.167*04 

.128*02 

.796-00 

.202-00 

.767*01 

-.675*01 

-.118*01 

•811*01 

-.253*01 

-.423*01 

9 

3705.44 

.186*04 

.840*01 

.437-00 

.312-00 

•S57*01 

-.457.01 

-.754-00 

.575*01 

-.167*01 

-.285*01 

10 

3622.98 

.201*04 

.555*01 

.203-00 

.369-00 

.356*01 

-.276*01 

-.428-00 

.361*01 

-.978-00 

-.171*01 

li 

3432.79 

.218*04 

.807*01 

.265-01 

.862-00 

.695-00 

-.332-00 

-.171-01 

.619-00 

-.833-01 

-.195-00 

12 

3346.96 

.219*04 

.328*01 

.595-05 

•100*01 

.230-03 

.565-04 

•497-04 

•135-03 

•594-04 

•461-04 

13 

3343.00 

.219*04 

♦204*01 

.887-06 

.100*01 

.238-06 

-.834-07 

.358-06 

-.477-06 

-.238-06 

-.917-08 

14 

3339.87 

.219*04 

.945-00 

.212-05 

.100*01 

.119-05 

.156-06 

.358-06 

•477-06 

.238-06 

•185-06 

15 

3335.78 

.219*04 

.263-00 

.652-06 

.100*01 

.143-05 

.218-06 

.477-06 

.715-06 

.119-06 

•401-06 

16 

3333.63 

•218*04 

.262-01 

.142-05 

.100*01 

-.715-06 

-.790-07 

-.238-06 

-.715-06 

-.119-06 

-.209-06 

17 

3333.36 

.218*04 

.288-03 

.256-06 

•100*01 

•  000 

-.109-06 

-.238-06 

•238-06 

-.119-06 

-.114-06 

-1 

3333.36 

.218*04 

.131-05 

.725-06 

•100*01 

-.238-06 

-.183-06 

-.119-06 

-.477-06 

-.238-06 

-.515-07 

Sample  Problem  Number  2 
Sample  Output 
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AEROThERM  EQUILIBRIUM  SURFACE  THERMOCHEMISTRY  JOB  2<IHI 

DERIVATIVE  PROPERTY  OUTPUT 

CP-FROZEN  CP-EOUIL  OLNM/ULNT  OLNM/ULNP  GAMMA 

.40311-00  .6ST82.05  -.38964-01  .16B3Z-0Z  .99838-00 

PROPERTY  ROUTINE  OUTPUT  IN  LB-MASS. F T«S£C .BTU. AND  OEG-R 

MU  1  MU 2  MOl.rT  HI1L  CPTIL  HTIL* 

•85S86-00  ,19801.02  .32790.0 2  .44.78.03  .70115-00  .61272.02 

ELEMENTAL  K  ANO  Z  MASS  FRACTIONS  BY  ATOMIC  NUMBER  ....  4 GAME*  a  .667) 
1  6  7  8  13  17 

.34892*01  .85266-01  .83931-01  .36347-00  .20000-00  .21245-00 

.89480-01  .13246-00  .13184-00  .33675-00  .39238-03  .30908-00 

SOLUTION  TIMES 

ITERATIONS  «  11  MACHINE  TIME  *  .411  SEC. 


CLOSED  SYSTEM  EQUILIBRIUM  SOLUTION  OUTPUT 

RELATIVE  MASaES  OF  COMPONENTS  1.2,  AND  3.  .70400-00  .96000-01  .20000-00 

MASS  CONOENSED/MASS  GAS  a  .60593-00 

THERMOOYNAMIC  state 

PRESSURE  <ATMOSPHERFS>  *  2.0000  TEMPERATURE  (DEG  K)  m  2318.000 

(LBF/IN2-ABS)  *  29.3920  1DEG  R)  *  4172.400 

GAS  CONDENSED  COMPOSITE 

DENSITY  (L8M/FT3)  =  .1340-01  .2152-01 

(GM/CM3)  a  .2147-03  .3447-03 

EMHALPY  (BTU/LEM)  •  3822*03  -.5901*04  -.2464*04 

(CAL/GM)  »  -.2123*03  -.3278*04  -.1369*04 

ENTROPY  CBTU/L8  OEG-R)  s  .2923*01  .7057-00  .2067*01 

MOLECULAR  WEIGHT  *  .2042*02  .1019*03  .3279*02 

CHEMICAL  STATE 

MOLE  FR.  *  MOLECULES  /  TOTAL  GAS  PHASE  MOLECULES 


SPECIES 

MOLE  FR. 

SPECIES 

MOLE  FR. 

SPECIES 

MOLE  FR. 

H 

.40973-02 

CO 

•21165-00 

N 

•10436-07 

0 

.34960-05 

AL 

.30759-07 

CLH 

.19360-00 

ALC 

.26533-21 

ALCL 

.34785-04 

ALCL2 

•39536-03 

AL20 

.38314-09 

C 

•16188-12 

CH4 

.99444-09 

C02 

.21122-01 

C2H2 

.14737-11 

C3 

.22070-21 

CL 

.19541-02 

CL2 

.16914-05 

HO 

.39032-03 

H2 

.30096-00 

H20 

.16746-00 

NO 

.11954-04 

N2 

.98226-01 

02 

•81884-06 

AL203* 

•12136*00 

C* 

•00000 

AL203* 

IS  .25452 

LIQUIO  AND  . 

74548  SOLID  8Y 

MASS. 

Sample  Preble*  Number  2 
Sample  Output#  Concluded 
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0120000000  2000.  1.0  1.0  AIR  PWP 

004 

7  NITROGEN  14.008  0*701 1 

fl  OXYGEN  16.000  0.2095 

16  ARGON  39.948  0.0094 

99  ELECTRON  0.00055 

1700000000000  OJANAF  03/61  N 

112965*6  134370*5  486'>44*1  383516-4  958460*5  480900*2  500.  3000.1  O.N 

112965*6  134370v5  428957*1  240844-3-417273*6  480900*2  3000.  5000.1  O.N 

1800000000000  OJANAF  06/62  0 

595590*5  135220*5  497228*1  380768-5  154749*5  500960*2  500.  3000.1  0.0 

595590*5  135220*5  657489*1-224268-3-891782*7  500960*2  3000.  5000.1  0.0 

1  18  00000000000  OCONVAIR  2PH-122  12/61  A 

000000*0*149400*5*418099* 1*122960-4*371200*7*484619*22000.  10000.1  A 

000000*0*149400*5*418099* 1 *1 22960-4* 371 200 *7*4846 19* 22000.  10000.1  A 

2700000000000  OJANAF  03/61  N2 

000000-0  221650*5  862699*1  116090-3-103715*7  637650*2  S00.  3000.1  0.N2 

000000-0  221650*5  984175*1-116232-3-612728*7  637650*2  3000.  5000.1  0.N2 

2800000000000  OJANAF  03/61  02 

000000-0  234460*5  804370*1  510872-3-152716*6  679730*2  500.  3000.1  0.02 

000000-0  234460*5  103071*2  290991-4-783079*7  679730*2  3000.  5000.1  0.02 

1718000000000  OJANAF  06/63  NO 

215800*5  227000*5  877623*1  899031-4-789656*6  688490*2  SOO.  3000.1  O.NO 

215800*5  227000*5  916260*1  657885-5-212519*7  688490*2  3000.  5000.1  O.NO 

2718000000000  OJANAF  12/60  N20 

195000*5  365450*5  144686*2  120489-3-153478*7  815900*2  500.  3000.1  0.N2O 

195000*5  365450*5  123036*2  459018-3  942382*7  815900*2  3000.  5000.1  0.N2O 

1  99  00000000000  OCONVAIR  ZPH-122  12/61  E- 

*149010*5*498851*1-272800-5-135900*6*164558*22000.  10000.1  E- 

♦149010*5*498851 ♦ 1-272800-5-1 35900*6* 164558*22000.  10000.1  E- 

1  18-1  99  000000000  OCONVAIR  2PH-122  12/61  A* 

363437*6*157330*5*486079* 1*21 1800-4* 165300*7*519257*22000.  10000.1  A* 

363437*6* 157330*5*486079*1 *21 1 800-4* l 6530 0*7*5 19257* 22000.  10000.1  A* 

17  -1  99  000000000  OCONVAIR  ZPH-122  12/61  N* 

*446641*6*151 310*5*50 1751*1 *61 7 100?4- 1841 00*  7, *696847* 220 06.  10000.1  N* 
*446641 *6* 151 31 0*5*50 175i *1*61 ^100-4- 1841 00* 7*496847*22000.  10000.1  N* 

1  7  r  8  -1  M  1  0  0  0  0  0  0  OCONVAIR  ZPH-122  I2'61  NO* 

♦23291 9*6*241 970*5*91 021 6*1 *277400-4-3 16600 *7 *654379*22000.  10000.1  NO* 

♦232919*6*241 970*5*9 l 02 16*1 *277400-4-3 16600 *7 *654379*22000.  10000.1  NO* 

18-1  99  000000000  OCONVAIR  ZPH-122  12/61  0* 

♦371 999*6* 149290 *5*336271 *1 *30671 0-3*590200* 7*484849* 22000.  1000C.1  0* 

♦37 1999*6* 149290*5* 336271 *1 *30671 0-3*590200 *7*484849*22000 •  10000.1  0* 


4000. 

6000. 

8000. 

10000. 


AIR  PROP 
AIR  PROP 
AIR  PROP 
AIR  PROP 


Sample  Problem  Number  3 
Listing  of  Input 
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AEROTMERM  EQUILIBRIUM  SURFACE  THERMOCHEMISTRY  JOB  AIR  PROP  OIZOOOOOCO 


RELATIVE  ELEMENTAL  COMPOSITIONS*  ATOMIC  WTS/UNIT 

MASS 

AT. NO.  ELEMENT  ATOMIC  WT 

COMPONENT  1 

;  COMPONENT  2 

COMPONENT  3 

7  NITROGEN  U.OQeQQ 

.0*32478 

-.0000000 

-.0000000 

8  OXYGEN*  16.00000 

•0140817 

-.0000000 

-.0000000 

18  ARGON  39.94800 

•0006400 

-.0000000 

-.0000000 

99  ELECTRON  ,00055 

•ooooooo 

-.0000000 

*.0000000 

I. 

70.  00. 

00.  00.  00*  00* 

OJANAE  03/61 

N 

•11296*06 

•13437*05  .48694*01 

a.  30352-04 

•95846*05 

•48090*02  3000.0000 

1 

.11296*06 

•13437*05  .42896*01 

.24084-03 

-.41727*06 

•46090*02  5000.0000 

1  .8016 

1* 

60.  00. 

0  0*  00.  00*  00* 

OJANAF  06/62 

0 

.59559*05 

•13522*0S  .49723*01 

.30077-OS 

.15475*05 

•50096*02  3000.0000 

1 

.59559*05 

.13522*05  .65749*01 

-.22427-03 

-.89178*07 

•50096*02  5000.0000 

1  *•*»» 

1. 

IB  0.  0  0, 

00.  00.  00.  00. 

OCONVA1R 

2PM-122  12/61 

A 

.00000 

.14940*05  .41810*01 

.12296-04 

.37120*07 

•48462*0210000.0000 

1 

.00000 

•14940*05  .41810*01 

.12296-04 

.37120*07 

•48462*0210000.0000 

1  1,25*3 

2* 

70.  00. 

00.  00.  00.  00. 

OJANAE  03/61 

N2 

•ooooo 

.22165*05  .86270*01 

.11609-03 

-•  10372.07 

•63765*02  3000.0000 

I 

.00000 

•22165*05  .98417*01 

-.11623-03 

-.61273.07 

•63765*02  5000.0000 

1  1.0807 

2. 

80.  00, 

00.  00.  00.  00. 

OJANAE  03/61 

02 

.00000 

.23446*05  .80437*01 

.51087.03 

-.15272*06 

•67973*02  3090.0000 

1 

.00000 

.23446*05  .10307*02 

•29099-04 

-.78300*07 

•67973*02  5000.0000 

l  1.1*66 

1. 

71.  8  0 . 

00.  00.  00.  00. 

OJANAE  06/63 

NO 

•21560*05 

.22700*05  .87762*01 

•89903-04 

-.78966*06 

•66849*02  3000.0000 

1 

•21580*05 

.22700*05  .91626*01 

.65788-05 

-.21252*07 

.68649*02  5000.0000 

1  132 

2. 

71.  80. 

00.  00.  00.  00. 

OJANAE  12/60 

N20 

.19500*05 

.36545*05  .14469*02 

.12049-03 

-.15348*07 

•61590*02  3000.0000 

1 

.19500*05 

.36545*05  .12304*02 

.45902-03 

.94238*07 

•81590*02  5000.0000 

1  i  <130 

1. 

99  0,  00. 

00.  00.  00.  00. 

OCONVAIR 

ZPH-122  12/61 

E- 

••00000 

.14901*05  .49685*01 

-.27280-05 

-.13590*06 

.16^56*0210000.0000 

I 

-•ooooo 

. 14901*05  .49885*01 

-.27280-05 

-.13590*06 

.16456*0210000.0000 

|*  ***** 

ELEMENT 

NITROGEN  OXYGEN 

ARGON 

ELECTRON 

BASE  SP 

N  0 

A 

E- 

1. 

18-1.  99  0. 

00.  00.  00.  00. 

OCONVAIR 

2PH-122  12/61 

A* 

•36344*06 

.15733*05  .40608*01 

.21180-04 

•16530*07 

.51926*0210000.0000 

1 

•36344*06 

.15733*05  .48608*01 

•21180-04 

•16530*07 

•Si 926* 021 0000 .0000 

1  1.25*2 

I. 

7-1.  99  0. 

00.  00.  00.  00. 

OCONVAIR 

ZPH-122  12/61 

N* 

•44664*06 

.15131*05  a 50175*01 

.61710-04 

-.18410*07 

.49685*0210000.0000 

1 

•44664*06 

. 1S1 31*05  .50175*01 

.61710-04 

-.18410*07 

•49665*0210000.0000 

1  .8016 

u 

7  1.  8-1. 

99  0.  00.  00.  0  0. 

OCONVAIR 

ZPH-122  12/61 

NO* 

.23292*06 

.24197*05  .91022*01 

.27740-04 

-.31660*07 

.65438*0210000.0000 

1 

.23292*06 

.24197*05  .91022*01 

.27740-0', 

-.31660*07 

.65438*0210000.0000 

1  1.1132 

I. 

8-1.  99  0, 

00.  00.  00.  00. 

OCONVAIR 

ZPH-122  12/61 

0* 

•37206*06 

•14929*05  .33627*01 

•30671- 03 

•59620*07 

.48465* 02 10 OOP «090Q 

1 

.37260*66 

.14929*08  .33627*01 

•30l67>J»3l 

.54020*07 

.4440S.0210000.0000 

1  .MO* 

aerotherm  equilibrium  surface  thermochemistry  solution  air  prop 


DERIVATIVE  PROPERTY  OUTPUT 

CP-FROZEN  CP-EQUIL  DLNM/DLNT  DLNM/OLNP  GAMMA 
,Z*6T0-00  .31T31-00  -.22871-02  .7AB3A-04  .12761*01 

PROPERTY  ROUTINE  OUTPUT  IN  LB-MASS.FT.SEC.BTU.ANO  DEG-R 

MU1  MU  2  MOL. XT  HTIL  CPTIL 

,10972*01  .26653*02  .29061*02  .BS0S6*03  .29751-00 


HTIL* 

.84985*03 


ELEMENTAL  K  ANO  Z  MASS  FRACTIONS  BY  ATOMIC  NUMBER 
7  8  18  99 

.74589-00  .22851-00  .25599-01  .00000 

.75412-00  .22250-00  .23375-01  .00000 


IGAMEX  •  .667) 


SOLUTION  TIMES 

ITERATIONS  •  22 


MACHINE  TIME  *  .677  SEC. 


CLOSED  SYSTEM  EQUILIBRIUM  SOLUTION  OUTPUT 

RELATIVE  MASSES  OF  COMPONENTS  1.2.  AND  3.  .10000*01  -.00000  -.00000 


THERMODYNAMIC  STATE 

PRESSURE  (ATMOSPHERES)  *  1.0000 

(LBF/1N2-A8S)  *  14.6960 


TEMPERATURE  (DEG  X)  *  2000.000 

(DEG  R)  *  3600.000 


DENSITY  (GM/CM3) 
(LBM/FT3I 


• 1T71-03  ENTHALPY  (CAL/GM)  *  471.348 

.1105-01  (BTU/LBM)  *  848.426 


MOLECULAR  VEIGHT 


29.0606  ENTROPY  (BTU/LBM-R)  »  .2136*01 


CHEMICAL  STATE 


a' 

MOLE  ER. 

*  MOLECULES  / 

TOTAL  GAS 

PHASE  MOLECULES 

SPECIES 

MOLE  ER. 

SPECIES 

HOLE  ER. 

SPECIES 

MOLE  FR. 

N 

,78499-09 

0 

.29907-03 

A 

.16622-01 

£- 

.95365-12 

N2 

.76976-00 

02 

•20341-00 

NO 

.79086-02 

N20 

.40650-06 

A* 

.20913-26 

N* 

.47971-31 

NO* 

.95365-12 

0* 

•82434-24 

<9 


Sample  Problem  Number  3 
Sample  Output 
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imooooos  1.0  C»  IN  *  JR 

!  6 
0. 

?.0  1.0  0,6  0.?  0.18  0.17 

003 

6  CARBON  12.011  -1.0 

7  NITROGEN  U.008  -0.767 

8  OXYGEN  16.000  -0.233 


3600000000000  OJANAF  12/60  C3 

189670*6  366220*5  146441*2  6?2536-4-1682?7*7  798410*2  500.  3000.1  0.C3 

189670*6  366220*5  1447R2*2  792232-4-646877*6  798410*2  3000.  5000.1  0.C3 

1618000000000  OJANAF  03/61  CO 

-264170*5  223570*5  865040*1  117021-3-898211*6  653700*2  500.  3000.1  O.CO 

-264170*5  223570*5  115496*2-424139-3-131563*8  653700*2  3000.  5000.1  O.CO 

2700000000000  OJANAF  03/61  N2 

000000-0  221650*5  862699*1  116090-3-103715*7  637650*2  500.  3000.1  0.N2 

000000-0  221650*5  984175*1-116232-3-612,128*7  637650*2  3000*  5000*1  0.N2 

1  6  0  0  0  0.  0  0  0  0  0-  0  OJANAF  03/61  *  C 

170886*6  135500*5  444A’33*1  228125-3  409830*6  492870*2  500.  3000.1  O.C 

170886*6  135500*5  412212*1  261908-3  262886*7  492870*2  3000.  5000.1  O.C 

2600000000000  OJANAF  9/61  C2 

198999*6  246990*5  776612*1  696081-3  185649*6  68551«*2  500.  3000.1  0.C2 

198999*6  246990*5  104162*2  566841-4-640205*7  685519*2  3000.  5000.1  0.C2 

4600000000000  OJANAF  12/60  C4 

242321*6  511230*5  205903*2  623436-4-257703*7  986760*2  500.  3000.1  0.C4 

242321*6  511230*5  210714*2-434895-4-404939*7  986760*2  3000.  SOOO.l  0.C4 

5600000000000  OJANAF  12/60  C5 

242374*6  656230*5  264706*2  806528-4-337125*7  111641*3  500.  3000.1  0.C5 

242374*6  656230*5  271156*2-580271-4-543183*7  111641*3  3000.  5003.1  0.C5 

1800000000000  OJANAF  06/62  0 

595590*5  135220*5  497228*1  380763-5  154749*5  5009*0*2  500.  3000.1  0.0 

595590*5  135220*5  65748(9*1-224263^3-091782*7  900960*2  3000.  SOOO.l  0.0 
280  0  080000000  OJANAF  03/61  02 

000000-0  234460*5  804370*1  510972-3-152716*6  679730*2  500.  3000.1  0.02 

000000-0  234460*5  103071*2  290991-4-783079*7  679730*2  3000.  5000.1  0.02 

30</OOOOOOOOOO  OJANAF  6/61  03 

341000*5  360230*5  135311*2  238739-3-60651 1 *6  860919*2  500.  3000.1  0.03 

341000*5  360230*5  139121*2  125613-3-980654*6  860919*2  3000.  5000.1  0.03 

1700000000000  OJANAF  03/61  N 

1 12965*6  134370*5  406944*1  303516-4  950460*5  480900*2  500.  3000.1  O.N 

112965*6  134370*5  428957*1  240944-3-417273*6  460900*2  3000.  SOOO.l  O.N 

1626300000000  OJANAF  03/61  C02 

-940540*5  365350*5  144559*2  210386-3-182392*7  798480*2  500.  3000.1  0.CO2 

-940540*5  765350*5  156451*2-391561 -4?60276fi*7  798480*2  3000.  5000.1  0.CO2 

3  6  2  8  0.  0  0  0  0  0  .♦  0'  8  OJANAF  12/60  C302 

-830000*4  652128*5  263775*2  103904-3-350368*7  112723*3  500.  3000.1  0.C3O2 

-830000*4  652120*5  263543*2  882201-4-287091 *7  112723*3  3000.  5000.1  0.C3O2 

1617000000000  OJANAF  12/66  CN 

111000*6  2 32498*5  655906*1  115326-2  479517*6  669760*2  500.  3000.1  O.CN 

111000*6  232490*5  988013*1  313855-3-649453*7  669760*2  3000.  5000.1  O.CN 

2627000000000  OJANAF  3/61  C2N2 

738699*5  511070*5  188740*2  559856-3-896873*6  985479*2  500.  3000.1  0.C2N2 

738699*5  511070*5  208204*2  630229-5-346865*7  985479*2  3000.  5000.1  0.C2N2 

4627000000000  OJANAF  3/61  C4N2 

127500*6  801239*5  297490*2  856052-3-152255*7  133063*3  500.  3000.1  0.C4N2 

127500*6  801239*5  327212*2  952999-5-541654*7  133063*3  3000.  5000. \  0.C4N2 

1718000000000  OJViAF  06/63  NO 

215800*5  227000*5  877623*1  899031-4-789656*6  688490*2  500.  3000.1  O.NO 

215800*5  227000*5  916260*1  657085-5-212519*7  688490*2  3000.  5000.1  O.NO 

1728000000000  OJANAF  06/63  N02 

801100*4  345800*5  137819*2  315611-4-133765*7  848090*2  500.  3000.1  0.NO2 

801100*4  345800*5  130154*2  172586-3  175340*7  048890*2  3000.  5000.1  O.NO? 

2710000000000  OJANAF  12/60  N20 

195000*5  365450*5  144686*2  128409-3-153478*7  815900*2  500.  3000.1  0.N2O 

195000*5  365450*5  123036*2  459018-3  942382*7  015900*2  3000.  5000,1  0.N20 

1600000000000  OJANAF  03/61  C» 

000000-0  144120*5  506075*1  953976-4-766621*6  121290*2  500.  3000.2  100. C« 

000000-0  144120*5  405134*1  291605-3  307202*7  121290*2  3000.  S000.2  100. C» 

1001000009  10.0  c  IN  AIR 

1001000009  100.0  C  IN  AIR 

8 


Sample  Problem  Number  4 
Listing  of  Input 
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AFROTHERM  EQUILIBRIUM  SURFACE  THERMOCHEMISTRY  solution  c*  in  aim 


derivative  property  output 

CP-FROZEN  CP-EQUIL  OLNM/QLNT  DLNM/QLNP  GAMMA 
.38619-00  .19303*01  -.12929*01  .76738-01  .11104*01 

PROPERTY  ROUTINE  OUTPUT  IN  LB-MASS «FT »5EC*bTU* AND  OtG-R 

MU1  MUp  MOL • wT  HT1L  CHTIL  HTIL* 

.11120*01  .27018*0?  .30715*02  .92709*04  .38476-00  .92126*04 

ELEMENTAL  K  ANO  Z  MASS  FRACTIONS  BY  ATOMIC  NUMBER  t  *  •  •  (GAMEX  ■  .667) 

6  7  8 

.66667-00  .25567-00  .77667-01 

.65033-00  .26319-00  .80483-01 

SOLUTION  TIMES 

ITERATIONS  *  15  MACHINE  TIME  *  .998  SEC. 


OPEN  SYSTEM  OUTPUT  (SURFACE  EQUILIBRIUM* 

BPRIME  PYROLYSIS  GAS  «  .00000  BPR I ME  CHAR  *  .20000*01 

SURFACE  SPECIES  IS  C« 

THCRMOOYNAMIC  STATE  ADJACENT  TO  Trtfi  SURFACE 

PRESSURE  (ATMOSPHERES)  *  10.0000  TEMPERATURE  (DEG  K)  *  4458.747 


(LBF/IN2-ABS) 

s 

146.9600 

(DEG  R) 

»  8025. 745 

DENSITY 

(GM/CM3) 

u 

.8394-03  ENTHALPY 

(CAL/GM) 

«  5073. 8S4 

(L9M/FT3) 

at 

.5240-01 

(BTU/LBM) 

>  0132.938 

MOLECULAR  WEIGHT 

at 

30.7146  ENTROPY 

(BTU/L9M-R) 

*  .2458*01 

CHEMICAL  STATE  AOJACENT 

TO 

THE  SURFACE 

MOLE  FR. 

*  MOLECULES  / 

TOTAL  GAS  PHASE  MOLECULES 

SPECIES 

MOLE  FR. 

SPECIES  MOLE  FR. 

SPECIES 

MOLE  FR. 

C3 

.29431-00 

CO  .14909-00 

N2 

.22837-00 

c 

.59404-01 

C2  .12460*00 

C4 

.22153-01 

cs 

.32923-01 

0  .10743-05 

02 

.10733-11 

03 

.22862-22 

N  .11973-02 

C02 

.43431-07 

C302 

.20517-09 

CN  • / 3254-01 

C2N2 

.10539-01 

C4N2 

.41560-02 

NO  .19430-06 

N02 

.34433-15 

N20 

. 18584-10 

C*  .00000 

PUNCHED  CARO  OUTPUT  (80  COLUMNS)  .... 

(  10.0000  .00000  2.000004458.7474  .000  5073.854  5073. 8S4  1  C*  C*  IN  A|R 


Sample  Problem  Number  4 
Sample  Output 
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aerotherm  equilibrium  surface  THERMOCHEMISTRY  SOLUTION  c*  in  air 


DERIVATIVE  PROPERTY  OUTPUT 

CP-FROZEN  CP-EQUIL  OLNM/DLNT  OLNM/OLNP  GAMMA 
.31142-00  .31145-00  -.44537-05  .7h059-05  .12949*01 

PROPERTY  ROUTINE  OUTPUT  IN  Lb-MASS.FT *SECi BTU.AND  DtG-P 

MU1  MU?  MOL «  W  T  HTlL  CPTIL  HTIL* 

•10807*01  .25924*0?  .20015*02  .45021*03  .31142-00  .46023*03 

ELEMENTAL  K  ANO  Z  MASS  FRACTIONS  BY  ATOMIC  NUMBER  ....  (GAME*  =  .667) 

6  7  8 

.14888-00  .65281-00  .19831-00 

.14888-00  .65280-00  .19832-00 

SOLUTION  TIMES 

ITERATIONS  =  15  MACHINE  TIME  =*  .344  SEC. 


OPEN  SYSTEM  OUTPUT  (SURFACE  EQUILIBRIUM* 

BPRIME  PYROLYSIS  GAS  *  .00000  BPRIME  CHAR  *  .17492-00 

SURFACE  SPECIES  IS  C* 

THERMOOYNAMIC  STATE  ADJACENT  TO  THE  SURFACE 

PRESSURE  (ATMOSPHERES)  =  10.0000  TEMPERATURE  (DEG  K)  *  2330.186 


(LBF/IN2-ARS) 

= 

146.9600 

(DEG  R) 

*  4194.335 

DENSITY 

(GM/CM3) 

s 

.1465-02  ENTHALPY 

(CAL/GM) 

*  255.699 

(LBM/FT3) 

s 

.9145-01 

(BTU/LBM) 

*  460.258 

MOLECULAR  WEIGHT 

= 

28.0147  ENTROPY 

(BTU/LRM-R) 

*  .2099*01 

CHEMICAL  : 

STATE  ADJACENT 

TO 

THE  SURFACE 

MOLE  FR 

.  =  MOLECULES  / 

TOTAL  GAS  PHASE  MOLECULES 

SPECIES 

MOLE  FR. 

SPECIES  MOLE  FR. 

SPECIES 

MOLE  FR. 

C3 

.19978-08 

CO  .34721-00 

N2 

.65276-OC 

C 

.14669-08 

C2  .24225-09 

C4 

.65724-12 

C5 

.96775-12 

0  .19122-09 

02 

.10660-13 

03 

.62259-27 

N  .13712-07 

C02 

.82019-05 

C302 

.57861-11 

CN  .16625-05 

C2N2 

.14820-04 

C4N2 

.22952-07 

NO  .36033-08 

N02 

.25489-17 

N20 

.53180-1? 

C«  .00000 

PUNCHED  CARO  OUTPUT  (80  COLUMNS)  .... 

(  10.0000  .00000  • 1 74922330# 1860  .000  255.699  255. 6*9  1  C*  C*  IN  AIR 


Sample  Problem  Number  4 
Sample  Output,  Continued 
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Seaaple  Problea  Huaber  4 
Staple  Output,  Concluded 


5-17 

0110  51)0.  .1  1.0  0.  0.  CM*  INPUT 

5 

1  HYDROGEN  1,000  -.09)0 

0  CARSON’  18.011  -.5470  -.488 

7  NITROGEN  14.000  -.7650  -.0190 

8  OXYGEN  16.000  -.8350  -.3410  -.873 

14  51LIC0N  88.090  -.839 


1  100000000000  OJANAF  18/60  H 

581080*5  134830*5  460883*1  555469-4  173095*6  308680*8  500.  3000.1  O.H 

581080*5  134830*5  308758*1  315085-3  989744*7  388680*8  3000.  5000.1  O.H 

1600000000000  0  JANAF  03/61  C 

170886*6  135500*5  4444)3*1  888185-3  409830*6  498870*8  500.  3000.1  O.C 

170886*6  135500*5  418818*1  861908-3  868886*7  498870*8  JOOO.  5000.1  O.C 

8700000000000  0  JANAF  03/6'  N8 

000000*0  881650*5  868899*1  116090-3*103715*7  637650*,  300.  3000.1  0.N8 

000000-0  881650*5  984 1 75*1-1 16838-3*618788*7  637650*?  3000.  5000.1  0.N8 

16180000,  000  70  0  0  JANAF  03/61  CO 

-864170*5  88*70*5  86*4*1  117081-3-898811*6  653700*8  500.  3000.1  O.CO 

-864170*5  883570*5  115496*8-484139-3-131563*8  653700*8  3000.  SOOO.l  O.CO 

1  14  00000000000  0  JANAF  18/68  SI 

106000*6  140180*5  457930*1  301346-3  148859*6  519980*8  500.  3000.1  O.St 

106000*6  140180*5  457930*1  301346-3  148859*6  519986*8  506.  3000.1  O.SI 

106000*6  140180*5  478318*1  141019-3  859988*7  519980*8  3000.  5000.1  O.St 

3116000000000  0J4NAF  18/68  CM3 

319400*5  434190*5  168763*8  401085-3-461803*7  786040*8  500.  3000.1  0.CH3 

319400*5  434190*5  804899*8-108088-3-114698*8  786040*8  3000.  SOOO.l  0.CH3 

1600000000000  0  JANAF  06/68  0 

595590*5  135280*5  497826*1  380768-5  154749*5  500960*2  500.  3000.1  0.0 

595590*5  135880*5  65?469*l-82486S-3-891708*7  500960*8  3000.  5000.1  0.0 

2«000ap0006%g  0  JANAF  03/61  08 

000000-0  8344&0*?  8(4370*1  510678-3-158718*0  679730*8  500.  3000.1  0.08 

000000-0  234460*5  103071*8  890991-4-763079*7  679730*8  3000.  SOOO.l  0.02 

181  14  000000000  OJANAF  09/63  OSI 

-848000*5  833590*5  891444*1  545951-4-458091*6  697470*8  500.  3000.1  0.0S1 

-848000*5  833590*5  850773*1  138710-3  109918*7  697470*8  3000.  5000.1  O.OSI 

1116170080060  OJANAF  03/61  CHN 

312000*5  355930*5  137083*8  552243-3-886955*7  758620*8  S06.  3080.1  O.CHN 

318000*5  355930*5  178695*2-895858-3-183671 *8  758620*8  3600.  5000.1  O.CHN 

1617000008000  OJANAF  18/68  CN 

109000*6  832496*5  655966*1  115326-2  479517*6  669766*2  566.  3600.1  O.CN 

109000*6  838490*5  988613*1  3131.55-1-646453*7  669760*8  3600.  5000.1  O.CN 

81860  0000000  0J(J1AF  03/61  C2H8 

541900*5  488570-  1  '.89960*8  7*90  44-^.409639.7 .849698*2  500.  3000.1  0.C2H8 

541900*5  488570*.  8>:7T-r*8  389048-3-645297*7  849690*8  3000.  5000.1  0.C8H8 

8687000030300  OJANAF  3/61  C8N8 

738699*5  511070*5  168740*2  ;5><.56-3-896873*6  985479*2  500.  3000.1  0.C8N8 

738699*S  511070*5  808804*2  830229-5-346665*7  985479*2  3600.  SOOO.l  0.C8N8 

2100000000000  OJANAF  03/61  H8 

008000-0  212100*5  711963*1  621950-3-712694*6  484650*8  500.  3000.1  O.H? 

000000-0  812100*5  681794*1  589854-3  865106*7  484656*8  3606.  5000.1  O.H? 

1  186000000000  ODUFF  BAUER  6/61  C?H 

1 17395*6  3496?0*5  134810*2  469100-3-187509*7  781140*2  500.  3000.1  0.C2H 

117395*6  349620*5  148516*2  109402-3-503862*7  781160*8  3000.  5000.1  0.C2H 

8600000000000  OJANAF  9/61  C2 

198999*6  246990*5  776618*1  696081-3  185649*6  685519*8  500.  3000-1  0.C8 

198999*6  246990*5  104162*2  566841-4-640205*7  685519*2  3000.  5000.1  0.C2 

1  1  16000  0  00000  OJANAF  03/61  CH 

142006*6  221300*5  886079*1  302211-3-100184*7  616120*2  500.  3000.1  0.e« 

142006*6  221300*5  707091*1  463281-3  552860*7  616120*2  3000.  5000.1  O.CH 

1  I  1  14  000000000  OJANAF  12/60  HS1 

114000*6  230106*5  679959*1  194499-3-978882*6  659960*2  500.  3000.1  O.MSI 

114000*6  230100*5  995466*1-335662-4-581070*7  659960*8  3000.  SOOO.l  O.HSI 
3600000000000  OJANAF  12/60  C3 

189670*6  366280*5  166441*8  688536-4-1688?’*7  798410*8  500.  3000.1  0.C3 

189670*6  ’66880*5  144788*8  798838-4-646877*6  798410*8  3000.  5000.1  C.C3 

1  1  3600000000  0  ODUFF  BAUER  6/61  C3H 

187703*6  489620*5  194464*8  508379-3-311933*7  988820*8  500.  3000.1  0.C3H 

187703*6  489680*5  199588*8  245687-3-638514*6  986686*8  3080.  SOOO.l  0.C3H 

2  14  0000000  0.  00  ' 0  OJANAF  12/68  S18 

130900*6  843436*5  691154*1  108588-3-354376*5  753940*8  500.  3000.1  0.SI8 

130900*6  243430*5  106398*8-190701-3-845149*7  753940*8  3000.  5000.1  0.SI8 

16000  00000000  OJANAr  03/61  C* 

000000-0  144180*5  586075*1  953976-4-766681*6  181890*8  500.  3000.8  100. C* 

000000-0  144120*5  485134*1  291065-3  307208*7  181890*8  3000,  5000.8  100. C® 

161  14  000600000  OJANAF  18/68  CSI* 

-806400*5  334819*5  158771*2-888607-3-503687*7  896880*8  500.  8773.8  100. CSI* 
-756200*4  348054*5  ’S8560*8-219419-4-170083*7  348989*2  8773.  5000.3  100. CSI* 
86  1  14  000000000  OJANAF  18/68  08SI* 

-2!7S00*6  468950*5  176588*8  549817-3-836409*7  473990*8  500.  1883.8  100.08SI* 

-215948*6  507890*5  817500*8-193988-4-888763*6  499160*8  1883.  5000.3  100.08SI* 
1  14  00000600000  OJANAF  12/68  SI* 

000000-0  178470*5  658874*1  838885-4-388164*6  185450*8  500.  1685.8  100. SI* 

180980*5  165680*5  600789*1  834683-4  356236*6  254880*8  1685.  5000.3  100. SI* 
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